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ABSTRACT 

We present an analysis of the chemical and ionization conditions in a sample of 100 weak Mg n 
absorbers identified in the VLT/UVES archive of quasar spectra. In addition to Mg II, we present 
equivalent width and column density measurements of other low ionization species such as Mg I, 
Fe II, Al II, C II, Si II and also Al in. We find that the column densities of C n and Si n are 
strongly correlated with the column density of Mg II, with minimal scatter in the relationships. The 
column densities of Fe II exhibit an appreciable scatter when compared with the column density 
of Mg II, with some fraction of clouds having 7V(Fe n) ~ A^(Mg n), in which case the density is 
constrained to njj > 0.05 cm~ 3 . Other clouds in which AT(Fe Ii) << iV(Mg n) have much lower 
densities. From ionization models, we infer that the metallicity in a significant fraction of weak Mg II 
clouds is constrained to values of solar or higher, if they are sub-Lyman limit systems. Based on 
the observed constraints, we hypothesize that weak Mg II absorbers are predominantly tracing two 
different astrophysical processes/structures. A significant population of weak Mg II clouds, those in 
which AT(Fe ii) << N(Mg ii), identified at both low (z ~ 1) and high (z ~ 2) redshift, are likely 
to be tracing gas in the extended halos of galaxies, analogous to the Galactic high velocity clouds. 
These absorbers might correspond to a-enhanced interstellar gas expelled from star-forming galaxies, 
in correlated supernova events. The AT(Mg ii) and AT(Fe n)/iV(Mg ii) in such clouds are also closely 
comparable to those measured for the high velocity components in strong Mg n systems. An evolution 
is found in iV(Fe n)/AT(Mg ii) from z = 2.4 to z = 0.4, with an absence of weak Mg II clouds with 
A^(Fe ii) ~ A^(Mg n) at high-z. The N(Fe n) ~ A^(Mg n) clouds, which are prevalent at lower 
redshifts (z < 1.5), must be tracing Type la enriched gas in small, high metallicity pockets in dwarf 
galaxies, tidal debris, or other intergalactic structures. 

Subject headings: galaxies: evolution — halo — intergalactic medium — quasars: absorption lines. 



1. INTRODUCTION 

The H i gas directly associated with galaxies that 
intercept the line of sight to background quasars ap- 
pears as optically thick Lyman Limit systems in the 
quasar spectrum. The prominent metal lines associated 
with these intervening absorbers are typically observed 
to be kinematically broad (Av ~ 100 — 400 km s" 1 ), 
strong, and often satura ted (e.g., Stcidel fe Sargenllll992t 
IChurchill fe Vogtl l200fi) . Studying the properties of a 
large population of such strong Mg II absorbers is a tech- 
nique used for constraining the evolution of metals in the 
interstellar media, gaseous halos a nd coronae of galaxies 
over a large history of the universe (|Lanzetta et al. Ifl987t 
IChurchill et al. Ill996f ). Apparently distinct from these 
strong Mg II absorbers are the population of quasar ab- 
sorption line systems in which the low ionization metal 
lines are weak. These systems are separated from the 
strong ones based on the standard definition of the rest- 
frame equivalent width of Mg II A2796 A line being 
Wr(2796) < 0.3 A. This is not a firm criterion for divi- 
sion, but has been followed as a convention on the follow- 
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ing basis. The survey of lSteidel fe Sargent! (|1992h . which 
identified a large population of strong Mg II absorbers 
used a sample of intermediate resolution spectra (AA ~ 
5 A) which had an equivalent width threshold of ~ 0.3 A. 
Later surveys of higher sensitivity and spectral resolution 
found that the equivalent width distribution of Mg II sys- 
tems at z ~ 1 increases steeply for W r (2796) < 0.3 A 
such that ~ 67% of all Mg II abso rbers (down to 02 A) 
from that epoch are in fact weak (|Churchill et al. Ill909fc 
iNaravanan et al.~l l2007). It later became clear that such 
an empirical basis for the classification of Mg II systems 
into strong and weak does bear some physical significance 
in that the two classes might b e tracing two or more 
different populations of objects (|Churchill et al. 1 119991 : 
IRigbv et al. Il2002t ICharlton et al. Il2003h . 

The class of weak Mg II quasar absorption systems 
have several remarkable properties that are unique. To 
begin with, unlike the strong systems, the weak Mg II 
systems are optically thin in neutral hydrogen and pro- 
duce metal lines th at are narrow [b(Mg) ~ 4 km s _1 ] and 
often unsaturated (jChurchill et al. Ifl9 99l . If weak Mg II 
absorbers are sub-Lyman limit systems with 10 15 8 < 
N(H i) < 10 16 8 cm -2 , they would account for a sig- 
nificant fraction (> 25%) of the high colum n density 
regime of the Lya forest (|Rigbv et al. Il2002j ). burveys 
of quasar fields to identify host galaxies have not often 
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found weak Mg n systems at close impact parameters 
(physical distance, D < 30 Kpc) of luminous star form- 
ing galaxies (L> 0.05 L*} (IChurchill fe Le Brunl 119981; 
IChurchill etaTI 120051 : iMilutinovic et al. I l2006f) ~This is 
a surprising result particularly in light of the fact that, 
in a substantial number of weak systems, the metallic- 
ity of the low ionization gas where the Mg n absorp- 
tion arises is constrained to values greater than O.IZ©. 
In some cases the best constraints require metallici- 
ties that are between and 10Z W 



ICharlton et al. 1 120031: iMisawa et all 



Rigb v et al. 1 120021 : 



2007). 



Thus, even 

though they have H I column densities that are ~ 4 or- 
ders of magnitude smaller than DLAs, weak Mg n ab- 
sorbers are produced in gas clouds with metallicities that 
are 0.5 - 2 dex higher than the average metallicity of DLA 
absorbers. 

The astrophysical systems associated with weak Mg II 
absorbers have not been identified yet. Several pos- 
sibilities exist, which partly account for the observed 
statistical and physical properties of weak Mg n sys- 
tems. Ex amples include extragal actic high veloc- 
ity clouds dNaravanan et al. I 120071 ). dwarf galaxies 
(|Lvnch et al. II2006T). gas clouds expelled in super winds 
from dwarfs (e g.lZonak et al. 1 12004 iStocke et al~ll2004l : 
iKeenev et aU l2006f ) and/or massive starburst galaxies 
and metal enriched gas in intergalactic star clusters 
(Rigb v~e"t al. I l2002t ). Recently, a number of authors 
have reported the detection of several high metallicity 
(> 0.1 Z Q ) gas clouds that are residing in the inter- 
galac t ic medium, both in the local universe (|Aracil et al. I 
2006t iTripp et al. I I2006D and at hig h reds hift [z > 2, 



Simcoe et al. I (|2004D : iSchave et al. I (f2007f) ]. The ob- 



served column densities of Mg II, C II, Si n and Fe II in 
these gas clouds, and the metallicities inferred for them 
are compa rable to several weak Mg II absorbers stud- 
ied so far (iRigbv et al. 1 [20021: ILvnch fe Charlton] [20071 : 
IMisawa et al. Il2007l) . 

Photoionization models of specific weak Mg II sys- 
tems have shown that they possess a two phase struc- 
ture. The low ionization gas, traced by such species 
as Mg II, Fe n, C n, Si n, etc. has a gas density of 
njj > 0.1 cm~ 3 that is roughly 2-3 orders of magni- 
tude higher than the density of t he associated high ion- 
ization gas traced by C IV lines (ICharlton et al. I [20031 : 
ILvnch fe Charlton! 120071 IMisawa et al. II2007D . The col- 
umn density measured for weak Mg II lines is typically 
N(Mg n) ~ 10 12 - 10 13 cm~ 2 . For such relatively small 
values, the high number density of ions derived from 
the photoionization modeling constraints the thickness 
of the low ionization gas to ~ 10 pc. The weak Mg n 
population occupies a significant volume of the universe 
with a cross-section similar to th e absorption cross sec- 
tion of luminous (-£*) galaxies ([Churchill et al. 1 119991 : 
iNaravanan et al.~| [2007). Thus, given their small thick- 
ness, if their gas clouds had a spherical geometry then 
they would be a million times more numerous than lumi- 
nous galaxies at z ~ 1 in order to reproduce the observed 
cross-section on the sky. However, an analysis comparing 
the relative incidence of high and low ionization gas in a 
sample of weak Mg n and C iv systems at low redshift 
(z < 0.3) favors a filamentary or shee t-like configuration 
for th e absorber's physical geometry ([Milutinovic et al~l 
2006) , instead of millions of individual, spherical low ion- 



ization Mg II clouds of ~ 10 pc size, embedded in a ~ kpc 
higher ionization halo, traced by C IV. 

In a ddition J Lynch et al. I (|2006l ) and lNaravanan et al. I 
(2007) discovered an evolution in the redshift number 
density (dN/dz) of weak Mg n absorbers over the inter- 
val 0.4 < z < 2.4. The dN/dz was found to peak at 
z = 1.2, and subsequently decline (from a no-evolution 
trend,) towards higher redshift. The equivalent width 
distribution function was also found to be different be- 
tween low and high redshift. At z ~ 1, the equivalent 
width distribution of weak Mg II is significantly higher 
than an extrapolation of the exponential distribution for 
strong Mg II absorbers. In contrast, at z ~ 2, the equiv- 
alent width distribution of weak Mg II clouds is only 
slightly in excess of an extrapolation of the strong Mg n 
distribution. In the context of these observed changes 
in the absorber statistics between z <~ 1 and z ~ 2, it 
becomes important to investigate if the changes are re- 
flective of an underlying change in the absorbers' phys- 
ical or chemical properties. Such an investigation may 
yield valuable clues into the physical nature of the kind 
of astrophysical processes/structures that produce weak 
Mg II absorption at low and high redshifts. The large 
sample size considered here, provides the scope for such 
an analysis. 

In this paper we present constraints on the chemical 
abundances, metallicity and ionization conditions for a 
sample of 100 weak Mg n absorbers identified in spectra 
that were extracted from the VLT/UVES archive. We 
compare the observed line properties between the various 
chemical elements in order to constrain the metallicity, 
density and line-of-sight thickness of the low ionization 
gas. Previous studies have focused on individual weak 
Mg II systems, whereas here we have attempted to de- 
rive the range of properties for a large ensemble of weak 
Mg II systems, which have only recently been discovered 
(Narayanan et al. 2007). In § [2] and § 02 we explain 
the measurement of line parameters and the comparison 
between the prominent metal lines in these systems. In 
§ 21 we present the Cloudy photoionization constraints 
on the densities and metallicities of these absorbers. An 
observed trend in the Fe n to Mg II ratio with redshift is 
presented in § [5j and a comparison of weak Mg n clouds 
with the high velocity subsystems in a sample of strong 
Mg II systems in § [6] We conclude the paper with a 
summary of the significant results (in § [7]) and a detailed 
discussion (§ [5]) on the nature of the gaseous structures 
selected by weak Mg n absorption. 

2. THE SAMPLE OF WEAK Mg n ABSORBERS 

The sample of weak Mg n systems presented in this 
study were identified in 51 quasar spectra extracted 
from the ESO archive. The quasars were observed 
using the UVES hi gh resolution echelle spectrograph 
dDekker et al. l[2000l) on VLT at R ~ 45, 000 (FWHM = 
6.6 km s _1 ). The detailed information on each quasar 
observation, such as the exposure time, emission redshift 
of the quasar, wavelength coverage, pr ogram ID and PI of 
the ob servation are listed in Table 1 of lNaravanan et al. I 
(2007). The reduction and wavelength calibration of the 
spectra were carried out using the ESO-provided MI- 
DAS pipeline. When multiple exposures of the same 
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target were available, they were co-added to enhance 
the S/N ratio of the spectrum, after applying vacuum- 
heliocentric velocity corrections. The final combined 
spectrum was continuum fitted using 1RAF 3 , and sub- 
sequently normalized. T he detailed reducti o n pro cedure 
can be found in § 2.1 of iNaravanan et all ((2007). The 
redshift path length searched in each quasar spectrum 
for weak Mg n AA2796, 2803 lines excluded regions blue- 
ward of the Lya emission to avoid contamination from 
the Ly-a forest. 

Th e 100 Mg II systems that we analyze here are taken 
from [Narayan an" et al. I (|2007l ) which described a survey 
for weak Mg n systems at 0.4 < z < 2.4. In addi- 
tion to magnesium, we use the information from ions 
of four other elements, viz. iron, carbon, silicon and 
aluminum, in each system to estimate the chemical and 
ionization conditions in the gas. Specifically, the lines 
that we consider are the following; Mg I A 2853 A; 
Mg ii AA2796,2803 A; Fe n AA2383,2600 A; 
Al in AA1855, 1863 A; Al n A 1671 A; C n A 1335 A; 
and Si n A 1260 A 4 . The coverage of the individual lines 
vary depending on the redshift of the system and the 
wavelength coverage of the spectrum in which the line is 
identified. The system plots are presented in Figure [1] 
(Note: The full set of system plots will be available in 
the online version of the journal. Here we provide only 
a few as examples) . 

Most weak Mg n absorbers also have associated high 
ioniz ation gas in a separate phase, traced by C iv lines 
(e.g. IChurchill et al. M 1999). The relative incidence of 
high and low ionization phase is a useful constraint 
on the physical geometr y of the absorber population 
(|Milutinovic et al. II2006T ). The UVES spectra offer si- 
multaneous coverage of C iv and Mg n over the redshift 
interval 1 < z < 2.4. Within this interval, in almost 
all cases C IV is detected. However, the C iv absorp- 
tion profile is in many cases offset in velocity with Mg II, 
indicating the presence of a separate phase for the high 
ionization gas. The C iv AA1548, 1550 profiles are shown 
in the various system plots of Figure Q] (Note: The full 
set of system plots will be available in the online version 
of the journal. Here we provide only a few as examples) . 
In this paper, our focus is on determining the ionization 
conditions and metallicity in the low ionization gas, and 
hence we defer the detailed analysis of the high ionization 
C iv phase and its association with the low ionization gas 
to a forthcoming paper. 

2.1. Measurement of Equivalent Widths 

For each system within the redshift interval of 
0.4 < z < 2.4, besides Mg n AA2796,2803 A, only 
Mg I A2853 A and Fe II AA2383, 2600 lines have wave- 
length long enough to be in the regions of the spec- 
trum typically uncontaminated by H I lines in the forest. 
The other prominent metal lines that we have measured 
have rest-frame wavelengths, A < 2000 A. As a conse- 
quence they are susceptible to blending with Lya forest 

3 IRAF is distributed by the National Optical Astronomy Ob- 
servatories (NOAO), which are operated by AURA, Inc., under 
cooperative agreement with NSF 

4 The wavelengths are vacuum wavelengths rounded to the near- 
est natural number 



lines, particularly since the redshift of the intervening 
absorber is often much less than the emission redshift of 
the quasar. In instances where line blending is appar- 
ent, we quote an upper limit on the measurement of the 
rest- frame equivalent width. For doublet/multiplet lines 
such as Fe n and Al in we have measured the equivalent 
width of the stronger member of the doublet. We also 
quote a 3er upper limit when a line is not detected at the 
3a level. Table 1 lists the rest-frame equivalent width 
measured for the various lines in each system. 

2.2. Measurement of Column Densities 

Absorption lines were fit with a Voigt function to es- 
timate the column density. An initial model for the line 
profile was derived using th e automated profile fitter AU- 
TOVP (|Dave et al. 111997ft . The AUTO VP routine gen- 
erated its model profile by performing a Voigt profile de- 
composition of the absorption feature and subsequently 
minimizing the % 2 by adjusting the velocity (v), column 
density (N) and Doppler parameter (b) for all the compo- 
nents in the model. The output of AUTOVP was then 
refined using a maximum likelihood least square fitter, 
MINFIT, which returns a best-fit model with a mini- 
mum n umber of Voigt profile components based on an 
F-test ((Churchill et al. Il2003f h To retain a component, 
requires an improvement in the model fit at an 80% sig- 
nificance. MINFIT derives the model absorption profile 
after convolving with a Gaussian kernel of FWHM = 
6.6 km s _1 , corresponding to the UVES spectral resolu- 
tion of R = 45,000. The column density and Doppler 
parameter with their la errors are obtained for this final 
model. 

Voigt profile fits were applied to the following lines 
associated with each system : Mg II AA2796, 2803; 
Mg i A2853 A; Fe n AA2383, 2600; Al m AA1855, 1863; 
Al II A1671 A; C n A1335 A; and Si n A1260 A. In the 
situation where a line is not detected at the 3a level, we 
quote an upper limit on the column density determined 
from the 3tr limit on the equivalent width. Our sample 
consists of only relatively high-S/N spectra. The 3a lim- 
its are hence low, so that we can assume linear part of 
the curve-of-growth for estimating the corresponding up- 
per limit in column density. To get robust constraints on 
the fit parameters, we use both members of the doublet 
while fitting profiles for Mg n and Al in, and both of 
the strong members of the multiplet in the case of Fe II, 
viz. Fe II AA2383, 2600. Weaker members of the Fe n 
multiplet were rarely detected at the 3-a level. By simul- 
taneous fitting of members of the doublet/multiplet, it is 
possible to recover the true column density, even if the 
stronger me mber of the dou blet/multiplet is saturated 
(see Sec 4.4.2 IChurchilllll997l ). Thus for example, in the 
case of Mg n, by using both members of the doublet, it 
is possible to recover the true column density for values 
ofjVYMg n) up to 10 14 cm" 2 (see Figure 4.3 IChurchilll 
1997). For lines that are not doublets, Voigt profile fits 
are unique only when the lines are unsaturated. In our 
sample, this would a problem only for the strongest of 
C ii A1335 A and Si n A1260 A lines. 

Table 2 lists the line parameters (v, N and b) thus 
measured for the various lines in each system. As men- 
tioned earlier, the lines with rest-wavelength A < 2000 A 
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are often found within the region of the spectrum that 
is contaminated by the Lya forest. For Al III, we found 
that blending with H I lines of the forest could be iden- 
tified by comparing the profile shapes of the individual 
members of the doublets. For the rest of the transitions, 
their profiles were compared to Mg n to rule out pos- 
sible contamination. Figure Q] shows the line profiles of 
the various low ionization transitions and Al in associ- 
ated with each system in our sample. For each line, the 
positions of the individual clouds, determined from Voigt 
profile fitting, are labeled. 



3. RESULTS FROM MEASUREMENT OF METAL LINES 

3.1. The Population of Single and Multiple Clouds 

From comparing the frequency distribution of the num- 
ber of clo uds per syst e m be tween strong and weak 
absorbers, iRigbv et al. I (|2002D discovered that unlike 
strong absorbers, weak Mg n systems have a non- 
Poissonian frequency distribution. Approximately two- 
thirds of the weak systems in their sample of 30 at z ~ 1 
had absorption in a single cloud, isolated in redshift. 
The clouds were narrow (b ~ 4 km s _1 ) indicating a 
small temperature and velocity dispersion in the gas. 
These systems were consequently called single cloud weak 
Mg II absorbers referring to the low ionization gas in 
a single narrow component, unresolved at R = 45, 000 
(FWHM= 6.6 km s _1 ). The other set of weak absorbers 
were called multiple cloud weak Mg n systems as they 
had the low ionization absorption in multiple clouds that 
are resolved at R — 45, 000 and kinematically broad 
(Av > 30 km s _1 ) compared to single clouds. 

The incidence of the number of low ionization clouds 
in any given weak Mg n system is important for consid- 
ering the physical geometry of the absorbin g structure 
(jEllison et al. Il2004t iMiTutinovic et al. Il2006h . Figure [2] 
shows the distribution of the number of Voigt profile com- 
ponents per system in our sample. In nine systems 5 , 
we found the Mg II line profile to have a slight asym- 
metry, sometimes yielding two components in the Voigt 
profile model. We have classified these as single cloud 
systems, and plotted them in the iV c = 1 bin, since the 
low ionization gas is predominantly still in a single gas 
cloud with an internal velocity dispersion is less than 
6.6 km s^ 1 . A similar asym metry in single c loud line 
profiles was also noticed by IChurchill et al. I (|1999h in 
HIRES/Keck high resolution spectra. However, their for- 
mal fitting procedure, with the lower S/N data, did not 
statistically favor a two component fit. The occasional 
asymmetry in the line profile is likely due to a contribu- 
tion to the low ionization absorption from a slightly off- 
set higher ionization gas cloud. Photoionization models 
have succeeded in reproducing the observed asymmetry 
in the line profile using a single low ionization phase and 
separate high ionization phases (e.g., see the ionization 
models for z = 1.405367 and z = 1.796237 systems in 
iLvnch fc Charlton! (|2007h and the z = 1.755704 system 

5 z = 0.599512 in Q2217-2818, z = 1.091866 in Q0042-2930, z = 
1.153704 in Q1151+068, z = 1.330502 in Q1157+014, z = 1.395635 
in Q0011+0055, z = 1.405367 in Q2347-4342, z = 1.491972 in 
Q0551-3637, z = 1.755704 in Q2243-6031 and z = 1.796233 in 
Q2347-4342 



in iMisawa et al. I (|2007f) . We have therefore chosen to 
classify the above nine systems as single cloud systems. 

Taking this into account, in our larger sample of weak 
systems we find that the single cloud absorbers ac- 
count for 48% of the total population, which is much 
smaller than th eir observed frac t ion in the HIRE S sam- 
ple de scribed in lChurchill et al.l (|1999h and lRigbv et al. I 
(2002), but is consistent within la. Within the red- 
shift interval of 0.4 < z < 1.4, identical to the re dshift 
path length covered by the IChurchill et al. I (|1999h sam- 
ple, we find that only 45% (34/76) of the weak absorbers 
are single cloud systems, indicating that our results are 
not affected by any evolutionary effect in which a larger 
fraction of z > 1.4 systems have multiple clouds. This 
is further confirmed in Figure [3] where we illustrate the 
distribution of single and multiple cloud absorbers as a 
function of redshift. We find weak absorbers showing 
absorption in both single and multiple clouds at all red- 
shifts within 0.4 < z < 2.4. A preference is not evident 
for a certain type of weak absorber (i.e. single or multiple 
cloud) towards either low or high redshift. 

3.2. Equivalent Width o/Mg n 

Figure 2] shows the distribution of rest-frame equiva- 
lent width of Mg nA2796 as a function of system red- 
shift. The strength of the low ionization phase as traced 
by Mg II demonstrates considerable scatter within the 
interval 0.4 < z < 2.4. A Spearmann-Kendall test sup- 
ports the null hypothesis that the equivalent width is 
statistically uncorrelated (Spearman's p — 0.04) with the 
redshift of the absorber. In strong Mg n absorbers, the 
low ionization absorption is never confined to a single 
cloud. The line profiles are often kinematically complex 
with the absorption spread in several clouds separated in 
velocity. From the statistical analysis of 23 strong Mg II 
system s along 18 quasar lines of sight, IChurchill et al. I 
(|2003f l found an average of ~ 8 clouds per system, with 
the absorption profile of one system resolved into as many 
as 19 different components. In addition, in that sam- 
ple of strong Mg n systems, a very strong correlation 
was found between the number of clouds and the rest- 
frame equivalent width IV r (2796). In the bottom panel 
of Figure [2 we illustrate that such a strong correlation 
(> 9a) also exists for weak Mg n systems, where most of 
the weaker systems (W,.(2796) < 0.1 A) have absorption 
only in one or two clo u ds. B oth IChurchill et al. I (|1999h 
and iNaravanan et al. I (|2007l ) found that the equivalent 
width distribution of weak systems at z ~ 1 rises rapidly 
towards smaller equivalent widths. This observation is 
also reflected in the bottom panel of Figure [2] where 
we find 67% of our sample of weak absorbers to be at 
W r (2796) < 0.15 A. 

3.3. Comparison of Rest-Frame Equivalent Width 

In FigureEJ we present the rest-frame equivalent width 
of the various metal lines compared to the equivalent 
width of Mg II A2796 line for both single and multi- 
ple clouds. The difference in the number of data points 
in each plot is attributed to the spectral coverage for 
the various transitions. A Spearman-Kendall nonpara- 
metric correlation test shows that the rest-frame equiv- 



5 



alent width of Mg I, Fe II, C II, Al II, and Al in are all 
correlated with the rest equivalent width of Mg n at a 
greater than 98% confidence level. Owing to fewer data 
points, Si II exhibits a correlation with Mg II of lesser 
significance (> 3cr). The Spearman and Kendall correla- 
tion tests were carried out using the ASURV astrostatis- 
tics package which takes in to consideration measure- 
ments that are uppe r limits (jFeigelson, fe Nelson] 119851 : 
lLavallev et"aT~lll992ft . 

We note that the C II equivalent width has a strong 
linear relationship with Mg II, with little scatter. The 
Si II may have a similar relationship, but it is hard to 
demonstrate with the smaller number of data points. All 
other transitions, though we have shown their equivalent 
widths to be correlated with Mg II, show a much larger 
scatter in this relationship. 

For Fe II, there is more than a factor of ten spread 
in the ratio ^(2383)/!^ (2796) at 0.2 < IU r (2796) < 
0.3 A. At small W r (2796), many of the W r (2383) mea- 
surements are upper limits, but a spread of more than a 
factor of two can still be demonstrated at Wr(2796) ~ 
0.05 A. A similarly large spread is also found for the rela- 
tionships between Mg I and Mg II, Al II and Mg II, and 
Al in and Mg II. 

The large scatter in the observed ratios between the 
various transitions can be brought about by a number 
of factors, and can be exploited to diagnose the physi- 
cal conditions of the absorber. For a given strength of 
Mg II, the spread in the strength of the other transitions 
can be due to variations in abundance patterns or to dif- 
ferences in the density/ionization parameters of the gas 
clouds (addressed in § 14.21 14.61 and [S]). There is also 
the possibility that absorption from two different ions 
arises in separate phases, in which case the scatter be- 
tween their equivalent widths could be quite large. If 
we are to distinguish between these different factors, it 
is important not to use the observed equivalent widths 
since they average together the contributions from dif- 
ferent clouds. The physical conditions are better probed 
through comparison of cloud-by-cloud column densities. 

3.4. Comparison of Column Densities 

In Figure [SJ we compare the Mg II column density 
measured for each cloud with the corresponding column 
densities in Mg I, Fe II, Si II, C II, Al II and Al in. In 
multiple cloud systems, the comparison is between each 
component of Mg II and the corresponding component in 
the other transition. Non-detections at the 3cr level are 
given as upper limits. To test for likely dependence be- 
tween the measured quantities, we apply Spearman and 
Kendall's non-parametric correlation tests. We find that 
the column densities of all ionization species except for 
Si II are correlated with the column density of Mg II at 
greater than 7a significance. The statistical measure of 
the correlation is smaller for Si II (3a significance) be- 
cause of fewer data, 40% of which are censored points. 
The strongest correlation is observed between N(C n) 
andiV(Mg n) (rank correlation coefficient, p = 0.762), in 
spite of being limited by fewer data points. Such a strong 
positive correlation in the equivalent width and column 
density of C II with Mg II justifies the use of C II lines, 



in conjunction with other low ionization lines such as 
Si II, to select weak absorbers at redshifts z > 2.5 where 
it becomes more difficult to use Mg n AA2796, 2803 lines 
because of the redshifting of these lines into the near- 
infrared regime. 

Among the various ions, the column densities of Fe II, 
C II and Si n display the least scatter with the column 
density of Mg II. The correlation between these ions and 
Mg II can be formalized as: 



log N(Fe n) = (1.15 ± 0.14) log N(Mg n) - 2.36, 
[a = 0.37) 

log iV(C n) = (0.82 ± 0.10) log N(Mg n) + 3.42, 
[a = 0.19) 

log AT(Si n ) = (1.02 ± 0.25) log iV(Mg n) - 0.26, 
[a = 0.28) 



The best-fit slope, the y-intercept and the correspond- 
ing la uncertainties of these regression lines were cal- 
culated using the survival analysis package AS URV 
Rev 1.2 (jFeigelson. fe Nelsonlll985l: lLavallev et al. 111992ft 
which implements the methods presented in llsobe et al. I 
(1986). The a values are the standard deviation of the 
respective fits. 

In Figure El we also plot the Solar composition of Fe, 
C, Si and Al with respect to Mg, for reference. The 
abun dance of carbon [log ( C/Mg^ = 0.970] is taken 
from lAllende Prieto et all 1)200 lU2002ft . of silicon [log 
(Si/Mg) = -0.0 301, iron [log (Fe / Mgk = -0.069] and 
magnesium from iHolweger et al. I (12001), and for alu- 
mmum [log (Al/Mg) Q = -1.110] from lGrevesse fe Sauval I 
(1998). The observed ratio of column densities between 
the various ions and Mg II when compared with the re- 
spective solar abundance ratios, can indicate whether the 
ionization fractions of C II, Si II, Al II, Al in and Mg I 
are comparable with that of Mg n in the low ionization 
gas. We note that this is the case for the observed Si II 
to Mg II and Al II to Mg II column density ratios, which 
closely follow the respective solar abundance ratios. On 
the other hand, the observed C n to Mg II ratios are 
above the solar abundance ratio, while the Fe n to Mg n 
ratios are below. A number of factors such as differ- 
ences in ionization parameter, differences in the elemen- 
tal abundances and/or contributions from different gas 
phases can combine to produce these observed trends, 
which are discussed in the next section. 



In addition, differential depletion of elements onto dust 
can lead to deviations from Solar composition. The pres- 
ence of dust has not been directly measured in weak 
Mg II systems. However, it has been found that dust 
extinction is significant o nly in stronger Mg I I absorbers 
(W r ( MK n) > 1.5 A; iKhare et alJ 120051 : lYork et all 

m 



2006). For low column density absorbers such as the 
weak Mg II systems, interstellar dust may not be a sub- 
stantial component influencing mctallicity estimates de- 
rived from gas phase abundances. In addition, CLOUDY 
models incorporating varying amounts of dust levels find 
dust having a negli gible effect on the density of the ab- 
sorbing gas as well (iRigbv et al.ll2002ft . 
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4. CHEMICAL AND IONIZATION PROPERTIES OF WEAK 
Mg n ABSORBERS 

The physical conditions in the low ionization gas clouds 
are constrained using the standard photoi onization code 
Cloudy (ver.07.02.01. lFerland et aI71ll998h . The primary 
objective is to derive limits on the metallicity, density, 
and line-of-sight thickness for the gas phase where the 
bulk of the Mg n absorption arises. For this purpose, the 
observed column densities of the other prominent low and 
intermediate ionization transitions - namely Mg I, Fe II, 
Si II, Al II, C II, and Al in, and their ratios to Mg II - 
are used. 

The ionization fraction for a given element is con- 
trolled by the density in the gas cloud as well as by the 
strength of the incident ionizing radiation. Weak Mg n 
systems are not known to reside at small impact param- 
eters (d < 30 kpc) from lumin o us sta r-forming galaxies 
(L > 0.05L*, IChurchill et~aTl (|1999l )). Hence the ion- 
ization balance in them is likely dictated by the inten- 
sity of the extr agalactic backg r ound radiation (EBR). We 
choose the lHaardt fe Madaul (|1996D model for the EBR 
which incorporates ionizing photons from quasars and 
star-forming galaxies after propagation through a thick 
IGM. A 10% escape fraction from galaxies is used for 
ionizing photons with A < 912 A. 

To determine the overall properties for our sample of 
weak Mg II absorbers, we generate a grid of Cloudy 
models for a range of ionization parameters (—8.0 < 
log U < —1.0) and neutral hydrogen column densities 
(14.0 < log A(H i) < 19.0). The weak Mg n systems in 
our sample span the redshift range 0.4 < z < 2.4. There- 
fore we consider two separate Cloudy grids modeled us- 
ing the integrated ionizing photon density (hv > 1 Ryd) 
at z = 1 and z — 2. The difference of ~ 0.5 dex in 
the intensity of the extragalactic background radiation 
field between these two redshifts does not critically af- 
fect the output of the Cloudy models. Nonetheless, to 
have a more tenable comparison between the data and 
the models, we plot the z < 1.5 and z > 1.5 systems on 
the 2 = 1 and the z — 2 grids respectively. We adopt a 
solar abundance pattern for the Cloudy models, but dis- 
cuss the effects of abundance variations. In the following 
sections, we discuss the constraints that the various ions 
provide towards the chemical and ionization conditions 
in the absorbing gas. 

4.1. Constraints from Mg II 

In our sample, the Mg n column densities of the in- 
dividual clouds in weak Mg II absorbers fall within the 
range 10 110 < N(Mg n) < 10 13 3 cm -2 . Our search for 
weak Mg II systems along the 81 quasar lines of sight 
is 86% complete down to the equivalent width thresh- 
old of W r (279 6) = 0.02 A, correspond ing to A(Mg n) ~ 
lO 118 cm" 2 ([Narayanan et al. I [2007). The N(Mg n) is 
useful to place limits on the metallicity of the low ioniza- 
tion gas phase. Figure [7] presents how the column den- 
sities of the various ions change with respect to the ion- 
ization parameter, log U, for different values of N(H i) 
and metallicity. For a given metallicity and ionization 
parameter, i.e., a certain density, an increase in A(H i) 
would correspond to an increase in the size of the ab- 



sorber. Also, with increasing ionization parameter, the 
neutral fraction of hydrogen declines such that to con- 
verge on the same value of N(H i), the size of the ab- 
sorber has to further increase. It is evident from Fig- 
ure [7| that for N(H i) = 10 15 cm" 2 , at sub-solar metal- 
licity (e.g. O.IZq), the model column density of Mg n 
is inadequate to explain the observed iV(Mg n) even for 
the weakest Mg n lines in our sample. For a given log 
U, the column densities of the ionization stages of var- 
ious elements scale almost linearly with both N(H i) 
and metallicity. Thus, higher A(Mg n) can be recov- 
ered by raising either iV(H i) or metallicity. For exam- 
ple, at N(H i) = 10 15 cm -2 , by raising the metallicity 
by 1 dex (to Z Q ), we find that the ionization models 
reproduce the observed column densities in the weaker 
Mg II systems (A(Mg n) < 10 12 cm" 2 ) in our sam- 
ple. For the same N(H i) value, at 10Z Q , a substantial 
fraction of the range of observed 7V(Mg n) is covered 
by the ionization models, except for those systems with 
N(Mg n) > 10 13 cm" 2 . Alternatively, with a 1 dex in- 
crease in A(H i), the curves shift correspondingly such 
that systems with A(Mg n) < 10 12 cm -2 can be pro- 
duced in 0.1Z Q gas. However, for N(H i) > 10 17 cm" 2 , 
the low ionization gas cloud is an optically thick, Lyman- 
limit absorber (i.e. able to produce a break in the spec- 
trum of the background quasar at A = 912 A in the 
rest- frame of the absorber). 

It can be concluded that the column density of Mg n is 
a suitable parameter for estimating limits on the metal- 
licity of the absorber. Our sample of weak Mg n systems 
spans a range of 2 dex in Mg n column density. Assum- 
ing solar abundance pattern, the metallicity in many of 
their low ionization gas clouds is constrained to be at 
least 0.1 Zq if the gas is optically thin in neutral hy- 
drogen (7V(H i) < 10 17 cm" 2 , see § [HJ>- Moreover, the 
strongest Mg n lines (A(Mg II) > 10 13 cm" 2 ) among 
the weak systems require supersolar metallicity. 

4.2. Constraints from Fe II 

In our sample of weak absorbers, 32% (66/205) of 
Mg ii clouds have Fe n detected at the > 3a level, out 
of which 81% are firm detections (i.e. detections unaf- 
fected by blending with other absorption features). The 
column density ratio, N(Fe n)/N(Mg n), falls between 
0.02 and 4.0. The range of values for the ratio remains 
unchanged even when we exclude upper limit measure- 
ments. Among the clouds with Fe n detected, 13 are 
single cloud systems and the remaining 53 are part of 
multiple cloud systems. 

Constraints for metallicity, similar to the ones derived 
using observed iVfMg n), can also be derived based on 
N(Fe n). Figure [7] illustrates how the column density of 
Fe ii changes with ionization parameter for different val- 
ues of A(H i) and metallicity. At N(U i) < 10 16 cm" 2 
and Z= O.IZq, A(Fe ii) < 10 112 cm" 2 which is inad- 
equate to explain the observed column density in sys- 
tems with Fe II detected. By raising A(H i) by one dex, 
we find a corresponding increase in the column density 
of Fe Ii in the models, such that a column density of 
A(Fe Ii) < 10 12 2 cm" 2 is possible at sub-solar metallic- 
ity. This still does not account for some fraction (~ 10%) 
of the observed Fe II lines. With metallicity increased to 
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solar and super-solar values the models begin to produce 
enough Fe n to explain the full range of observed values. 
In general, we can infer that for the systems in which 
Fe II is detected in our sample, the metallicity is con- 
strained to values of Z > Z Q if N(U i) < 10 17 cm -2 . 

It is evident from the column density comparison 
in Figure [6] that, for a given iV(Mg n), the observed 
iV(Fe n) has a spread of ~ 1 dex between the various 
systems. This spread is also evident in Figure [5] which 
compares the rest-frame equivalent widths. For gas that 
is optically thin, the ratio of column density between 
various ionization stages does not depend on metallicity 
since all individual column densities scale linearly. An 
exception to this can occur (discussed in § I4.4[) for cer- 
tain ions at supersolar metallicities where cooling leads 
to much lower gas temperatures. In the optically thin 
regime, for a given abundance pattern, the ratio of col- 
umn densities between different elements varies primar- 
ily with ionization parameter. The relative strength of 
Fe II compared to Mg II in a system is particularly sen- 
sitive to ionization parameter for log U > —4.0. Thus 
we over-plot, in Figure [8j the observed column density 
ratios of Fe n to Mg II on a Cloudy grid of photoioniza- 
tion models. The Cloudy grid is for a range of log U and 
iV(H i) at sub-solar, solar and super-solar metallicities. 
The censored data points that occupy the left of Figured] 
are systems in which Mg n is very weak. The Fe II, be- 
ing even weaker, is not detected at the 3cr significance 
threshold. The S/N of the best of our sample of quasar 
spectra are comparable and therefore the envelope of the 
ratio iV(Fe ll)/iV(Mg n) for censored data points is seen 
as increasing with decreasing A^(Mg n). 

We can place constraints of log U assuming that the 
Fe II and Mg n arise in the same phase. To begin with, 
we notice that the column density ratios of all clouds 
in our sample confine the ionization parameter to log 
U < —2.0, corresponding to uh > 0.002 cm~ 3 (for log 
n 1 = —4.70 at z = 2). At 0.1Z Q , the systems with 
Fe II detected require H I column densities greater than 
10 16 cm~ 2 . By increasing the metallicity, the grids shift 
to the right proportionately such that the same Fe II to 
Mg II ratios can now be recovered from weaker H I lines 
(JV(H i) ~ 10 14 -10 15 cm" 2 ). Thus, if the low ionization 
gas is thin in neutral hydrogen, the metallicity in systems 
where Fe II is detected is constrained to solar or super- 
solar values. 

In the sample of 17 single cloud weak absorbers stud- 
ied by Rigby et al. (2002), a subset of systems with 
log [iV(Fe n)/JV(Mg n)] > -0.3 were classified as iron- 
rich. These systems were found to have high metallic- 
ity (> O.IZq), particularly strong constraints on den- 
sity (log U < —4.0, tih > 0.09 cm -3 ), and small sizes 
[iV(H i) + iV(H n) < 10 18 cm -2 , R < 10 pc]. The rel- 
atively high Fe n to Mg II ratio indicated that the high 
density, low ionization gas in the iron-ri ch systems is 
not q- enhanced. Following the definition of iRigbv et al. I 
pool , we find that 30 clouds in our sample are iron-rich, 
excluding censored data points. Comparing the data to 
the Cloudy grid, we find that the ionization parameter 
in these systems is constrained to an upper limit on the 
ionizing parameter between -3.2 and -3.7 depending on 
the difference in ionizing photon number density between 



z = 2 and z — 1 respectively. A limit of log U < —3.7 
translates to a density of nn > 0.05 cm -3 in the low ion- 
ization gas for log ri 7 = —5.04 cm -3 (the number density 
of ionizing photons with hv > 13.6 eV at z = 1). The 
density constraint for the Fe n rich systems translates 
into a small upper limit for the thickness (R < 10 pc) of 
the absorber. In systems where Fe n is weak compared 
to Mg II, the constraint on density is much lower (log 
U < -2.0). 

In this analysis, we have assumed a Solar abundance 
pattern. Changing the abundance of any element from 
this pattern would result in a corresponding change in 
all ionization stages of that element. Thus, the iron- 
rich systems can have a lower constraint on density if 
the abundance of iron in the low ionization cloud is en- 
hanced relative to the solar abundance pattern, since the 
Cloudy grids would be shifted upwards. Such an abun- 
dance pattern is not physically well motivated. On the 
other hand, an a— enhanced abundance pattern is ruled 
out for these iron-rich systems as its effect would be to 
shift the Cloudy grids further down such that the ioniza- 
tion models will not be able to reproduce the observed 
Fe II to Mg II ratio. The a— enhancement is, however, 
conceivable for the clouds in which Fe n is low compared 
to Mg II. The ionization model, in that case, would infer 
higher densities for the low ionization gas. 



Finally, IRigbv et al71 (|2002f) found that the 
N(Fe ll)/iV(Mg II) in their HIRES sample had a 
bimodal distribution with an apparent gap of ~ 0.5 dex 
at -0.8 < log [N(Fe n)/JV(Mg n)] < -0.3. It was 
therefore used as basis for defining the iron-rich sys- 
tems, and to suggest that there may be a separate 
class where Fe II is weak relative to Mg II. Figure [5] 
shows the histogram distribution of the Fe n to Mg n 
ratio for our sample of weak Mg n single and multiple 
clouds. The bin size is equiv a lent t o the gap in the 
distribution that IRigbv et al. I (|2002l ) found for their 
sample. The distribution from our sample does not 
suggest a bimodality, either for single or multiple clouds. 
This remains true for smaller bin sizes as well. Hence 
the ap parent gap that was suggested by the IRigbv et al. I 
(2002) data can be attributed to inadequate sample size. 
Additionally, we also note that there is no difference in 
the observed Fe n to Mg n ratio between single and 
multiple cloud systems. The individual clouds in the 
multiple cloud systems have the similar log U constraints 
as single clouds, with iron-rich systems detected in both 
category. 

4.3. Constraints from Mg I 

In this section, we explain the constraints that are 
available from the observed Mg I to Mg n ratio in 
weak systems. In the past, single p hase photoion- 
ization models (using CLOUDY 90; iFerland et al. I 
(1998)) have failed to reproduce the obse rved Mg I to 
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Mg I I ratio in some strong Mg II systems ( 
12001 iChurchill et all 12003 IDing et al. I 
Mg i/Mg II ratio derived from the models was lower than 
the observed neutral to singly-ionized ratio. To circum- 
vent this, a separate phase was proposed in which the 
Mg I ionization fr action is higher (|Churchill et al. 1120031 : 
IDing et aT~ll2003f) . This separate phase would have a 



higher density (njj > 1 cm -3 ) and lower temperature 
(T < 600 K) than the gas phase associated with the 
Mg II absorption. The Mg I lines corresponding to such 
low temperatures are very narrow (6 ~ 2 km s -1 ) and 
are therefore unresolved at the R — 45, 000 of the earlier 
HIRES and UVES observations. However, through su- 
perhigh spectral resolution observations, at R = 120, 000 
(Av = 2.5 km s -1 ), it has been demonstrated that the 
Mg I lines are not narrower than what is derived for 
R = 45, 000 (Narayanan et al. 2007). 

Compared to Fe II, only a few weak Mg n systems in 
our sample have Mg i detected at the 3<r level. Out of 
the 200 weak Mg II clouds for which there is coverage, 
Mg I is detected in only 7 single cloud systems and in 20 
clouds in multiple cloud systems. Both single and mul- 
tiple clouds span roughly the same range of values for 
the Mg I to Mg II ratio, between -2.2 and -0.5, consider- 
ing only firm detections. The neutral magnesium fraction 
(Mg°/Mg fota /) is thus small, compared to the Mg II frac- 
tion (Mg n/Mg tota /), in these systems. This most likely 
explains the large scatter in the range of limits, evident 
in Figures [5] and [5] , since we are sampling a large num- 
ber of quasar spectra with differences in sensitivity. The 
spectra with the highest S/N ratio in our sample, how- 
ever, constrain the Mg I column density to values as low 
as 10 9 - 5 cm -2 , ~ 2 dex smaller than N(Mg Ii), indicat- 
ing that the neutral fraction in the Mg II phase is indeed 
not very high. 

Figure [TU1 shows the observations compared to the grid 
of Cloudy (version 07.02.01) ionization models. To be- 
gin with, the ionization models are able to recover the 
observed Mg I to Mg II ratio from a single phase. Com- 
pared to the CLOUDY (version 90) grid of ionization 
models presented in Churchill et al. (2003), the mod- 
els displayed in Figure [T0l have the Mg i/Mg n fraction 
higher by ~ 0.5 dex for a given log U. The difference in 
the ionization fraction of magnesium is a result of im- 
provements in the rate coefficients for charge transfer 
reactions, i ncorporated into the more modern versions 
of Cloudy (Kingdon et al. 1996). The relevance of the 
charge transfer reaction (H + Mg + — > H + + Mg) in con- 
trolling the ionizat ion fractions of Mg I and Mg II has 
also been noted bv lTappe fe Blackl ((2004). For our sam- 
ple of weak Mg II systems, the ionization models derived 
from the revised version of the photionization code sug- 
gests that a single phase solution is possible. In fact, the 
observed ratio of Mg I to Mg n in strong Mg II absorbers 
can also now be explained without invoking a separate 
cold phase for Mg I. 

We find that, in a large majority of the systems for 
which information on Mg I is available, the ionization 
parameter is confined to log U < —2.5, for solar and 
super-solar metallicity. At Z < Zq, the constraint on 
ionization parameter is higher by ~ 1 dex. The frac- 
tion Mg i/Mg II is expected to decrease with an increase 
in the ionization conditions in the gas. This is evident 
in Figure [7] Therefore, the systems with higher Mg I 
to Mg II ratios will have lower constraints on ionization 
parameter (log U < —3.0) and correspondingly higher 
constraints on density (jih > 0.02 cm -3 ), identical to 
iron-rich systems. Moreover, if the H I lines are weaker 
(iV(H i) < 10 16 cm -2 ), solar or supersolar metallicities 



will be necessary to reproduce the observed Mg I and 
Mg II column densities (see Figure [101) • 

4.4. Constraints from C II 

In our sample, all the 15 weak Mg II systems for which 
there is coverage of C II A1335 A show prominent C II 
lines. Among these, eight are single cloud and the re- 
maining are multiple cloud systems. The ratio of N(C Ii) 
to N(Mg ii) is always significantly greater than 1, with 
a range of values between 4 and 60. A larger oscillator 
strength for the Mg nA 2796 line, and a longer wave- 
length compared to C nA 1335, leads to them having 
comparable equivalent widths. 

The Cloudy grid of single phase photoionization mod- 
els, comparing the column density of C II to Mg II, is 
shown in Figure [11] We find that for most of the sys- 
tems, the ionization parameter would be constrained to 
log U > —2.5, implying a gas phase density of Uh < 
0.006 cm -3 . Such large values of log U would be incon- 
sistent with that inferred for the low ionization phase 
for the clouds in which Fe II is detected. However, in 
Figure [TBI we compare the observed C II to Mg II ratio 
against Fe II to Mg II in systems with simultaneous cov- 
erage of both lines. Only for two systems [z — 1.585464, 
and z — 1.988656), do we have firm (i.e. measurements 
that are not limits) detections for both C II and Fe II. 
In these two systems we find 7V(Fe n) to be ~ 1 dex 
smaller than 7V(Mg n), and therefore a constraint of log 
U > —2.5 (see the z > 1.5 iron grid). This is consis- 
tent with the log U derived using the C n to Mg II ratio 
in the corresponding systems allowing for a single low 
ionization phase solution. The density of this low ion- 
ization phase would be smaller than what is estimated 
for the iron-rich systems. Unfortunately, such a definite 
statement cannot be extended for all the other systems 
plotted in Figure [T5l as their Fe II measurements are up- 
per limits. Nonetheless, even these limits are consistent 
with a low density. 

Many of the detected clouds in Figure [8] did not ap- 
pear in Figure [TBI because C II could not be measured for 
them due to contamination or lack of coverage. For these 
clouds the inferred ionization parameters range from log 
U ~ -2.2 to ~ -3.5. If the C II and Mg n arise in 
same phase, we would expect the C II in the clouds in 
which Fe II is detected to have N(C ll)/JV(Mg ii) > 10, 
such that consistent log U values would be derived from 
C II. Such a direct comparison between the measured 
C II and Fe II may not apply because C II may arise 
partly from a higher ionization phase. In addition to 
the low ionization phase, most weak Mg II systems also 
have an associated high ionization phase where the den- 
sity is low (njj < 10 -3 cm -3 ). Although dominated by 
higher ionization states of carbon (C in and C iv), the 
C II ionization fraction (i.e. C u/Ctotai) can be non- 
negligible in this phase. For example, at log U = —1.5, 
iV(H i) = 10 15 cm -2 , and Z = O.3Z (typical values de- 
rived from photoioni zation models, e.g. see Table 5 of 
iMisawa et al~l (12001 ). N(C „) = io 13 - 3 cm" 2 , which is 
comparable to the detected C II in many of our clouds. 
The iV(Mg ii) contribution from this high ionization 
phase is negligible. In summary, even though C II is de- 
tected in all weak Mg II systems, since it does not arise 
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exclusively in the low ionization phase, it may not pro- 
vide as robust a constraint on the ionization parameter 
as does Fe II. 

We note, also, that the grid with a metallicity of 10Z© 
does not cover many of the C n data points. At this 
metallicity, temperatures fall to T < 1000 K because of 
metal cooling, even at log U > —1.5. The gas, includ- 
ing both magnesium and carbon, is less heavily ionized 
at the low temperatures, but the effect is stronger for 
Mg II so that the density of Mg n is larger relative to 
C II. Clouds that have supersolar metallicity constraints, 
based on other transitions, would then need to have a 
large contribution to the C n from a separate phase. 

4.5. Constraints from Si II 

In our sample, we find that the Si II column density 
is comparable to the column density of Mg II, as shown 
in Figure [H The ratio of column densities has values 
between 0.2 and 3.2, with majority of them at ~ 1. Most 
of the clouds have Si II to Mg II ratios that fall on the 
grid of Cloudy models in Figure [T2j However, the grids 
do not provide much leverage in determining log U since 
Si II and Mg n are similar over most of the parameter 
space, particularly for solar and higher metallicities. 

The single phase ionization models suggest that sys- 
tems in which AT (Si n) < A(Mg n) require log U > —1.5 
and high metallicity, assuming a Solar abundance pat- 
tern. This is evident in the Cloudy grids where data 
points corresponding to low Si II to Mg n column den- 
sity ratio are below the model expectations for Z < Z Q . 
Only for supersolar metallicities, do any of the models re- 
produce the low Si II to Mg II ratio. This is also evident 
in Figure [7] where the Si n and Mg II ionization curves 
cross only at super-solar metallicity for log U > —1.5, 
corresponding to Hh < 10~ 3 cm" 3 . The low Si II to 
Mg II ratio can also be obtained from single phase models 
by lowering the abundance of silicon compared to other 
a— process elements, in which case the metallicity could 
be lower, and the density higher. 

In general, the Si II does not provide a robust con- 
straint on the ionization parameter for the low ionization 
gas. In the small number of clouds for which both Si II 
and C II are covered they usually provide consistent con- 
straints on log U , taking into account that some fraction 
of the C II can arise in a separate phase. 

4.6. Constraints from Al II 

Figure [13] shows the Cloudy grid of photoionization 
models, with measurements of Al II to Mg II overplotted. 
Within the sample of systems, we find that A(A1 n) is 
always smaller than A(Mg n), sometimes by as much as 
0.5 to 1.5 dex. For a solar or smaller metallicity, many 
systems (those with iV(Al n) > 0.1A(Mg Ii)) are cov- 
ered by the photoionization grids. Particularly at so- 
lar metallicity, the ratio of Al n to Mg II is not very 
sensitive to ionization parameter, so it cannot be used 
to effectively measure this quantity. Of special note 
are a number of systems, both at high and low red- 
shift, with A(A1 ii) < 0.1A(Mg n). These points are 
not covered by the grid for all metallicities, based on a 



solar abundance pattern. The IOZq grid does extend 
(for \ogU > —1.5) to somewhat lower values of Al II 
to Mg II. A supersolar metallicity could help to explain 
some systems with a low Al II to Mg II ratio, such as the 
z = 1.68079 system towards Q 0429 - 4091 described in 
Misawa et al. (2007). However, for many of these sys- 
tems, the most likely explanation would be a reduction 
of the aluminum abundance relative to magnesium by up 
to ~0.7 dex. A reduction would effectively shift the grids 
down by the same amount. Such an abundance pattern 
is feasible, since it is consistent with a-enhancement. We 
note that the same shift should apply for these systems 
in the Al in grids as well. 

4.7. Constraints from Al III 

Among the 74 systems with simultaneous coverage of 
Mg ii AA2796, 2803 A and Al in AA1855, 1863 A lines, 
Al in is detected in 12 single cloud systems and in 44 
clouds in 20 multiple cloud systems. The ratio of Al III 
to Mg II column density falls within the range 0.04 to 
0.98, excluding upper limits. The Cloudy grids of single 
phase models are shown in Figure ITU Assuming that the 
Al in and Mg II are produced in the same phase, for the 
Al III detections, the ionization parameter ranges from 
-3.5 < log U < -1.0. 

In the previous section, in order to reconcile the Al II 
to Mg II grid with the points below that grid, we pro- 
posed a reduction of the aluminum abundance relative 
to magnesium. If this is applied to the Al III, it shifts 
this grid downwards so that it does not cover some of the 
data points. If the systems with small Al n to Mg II also 
have small Al in to Mg II this is not a problem, though 
it does require large log U for even these systems. This is 
found to be the case in Figure [T5l where we have plotted 
the ratio of Al in to Mg II versus Al II to Mg II. For 
all of the systems below the Al II grid (see Figure H"3"|). 
A(A1 in) < 0.3A(Mg ii). 

Although we have not identified specific clouds for 
which Al in and Mg II cannot arise in the same phase, 
we note that some of our Al in detections imply large 
ionization parameters, log U > —2.0. This is even more 
the case if we rely on a decrease of the aluminum abun- 
dance. So either there is a sub-population of weak Mg II 
clouds which are of a higher ionization state, or some of 
the Al in is produced in a higher ionization phase, such 
as the one giving rise to the bulk of the C IV absorption. 

4.8. Al in to Al II ratio 

In damped Ly-a (DLA) systems, the chemical abun- 
dance estimations are often carried out under the as- 
sumption that the ionization corrections are not signif- 
icant, since the gas is expected to be predominantly in 
the low ionization phase. However, the detection of Al III 
lines at the same v elocity as the l ow ion ization lines in 
several DLAs lead IVladilo et aTl (|2001l ) to investigate 
the relevance of i onization corrections for these systems. 
In their analysis. IVladilo et al. I (|2001f ) observed that the 
iV(Al m)/iV(Al ii) ratio in DLA systems exhibits an 
anti-correlation with A(H i). This relationship was de- 
scribed as intrinsic to DLAs, and was used to suggest that 
the Al in to Al II ratio in these systems could be a sensi- 



10 



tive probe of the ionization co nditions in the gas. Using a 
sampl e of sub-DLA systems, iDessauges-Zavadskv et all 
( 2003) examined if this anti-correlation extends to lower 
N(K i). They found the Al in to Al II ratio in sub-DLAs 
to be in the same range as for DLA systems. In other 
words, the anti-correlation trend did not seem to extend 
to sub-DLAs (10 19 < N(R i) < IP 2 03 cm" 2 ). However, 
more recently IMeiring et al. I (|2007f ) found that the anti- 
correlation could apply even to sub-DLAs, based on a 
different sample of systems. 

In our sample, 28 weak Mg II clouds have measure- 
ments of both Al III and Al II, of which 15 are firm de- 
tections in both (i.e. measurements that are not lim- 
its) . In Figure 1161 we plot their ratio with respect to 
the corresponding A^(Mg n). We also plot, in an ad- 
jacent panel, the Al in to Al II ratio in DLA and sub- 
DLA systems as a function of N (H i), based on informa- 
tion extracted from the literatur e dyiadilo e t al. 120011 : 
IDessauges-Zavadskv et all 120031: IMeiring et al. I 120071 ). 
We find the Al in to Al n ratio in weak Mg ii systems 
to be considerably higher than in DLAs and sub-DLAs. 
On average, the ratio is ~ 0.5 — 1 dex higher than what 
has been measured for the other two classes of systems. 
This indicates that the ionization conditions are higher 
in weak Mg II systems than in DLA or sub-DLA systems. 

Based on photoionization modeling (see § 14. 1[) we have 
concluded that a large fraction of weak Mg II clouds have 
a metallicity of solar or higher if the clouds are optically 
thin in neutral hydrogen. The observed redshift number 
density of weak Mg n absorbers is too large for all of 
them to be Lyman limit systems (as explained in § [5]). 
If N(R i) < 10 17 cm~ 2 for the weak systems plotted 
in Figure [T751 then we conclude that the anti-correlation 
trend discovered by Vla dilo et al. (200 1 ) for DLA sys- 
tems, and supported by IMeiring et al. I (|2007l ) for sub- 
DLA systems, continues to lower A^(H i) values. 

5. EVOLUTION OF THE LOW IONIZATION PHASE 
STRUCTURE 

One of our objectives in carrying out the chemical and 
ionization analysis on a large sample of weak Mg II sys- 
tems is to find out if there are any evolutionary trends 
observable in the absorber population. Our VLT/UVES 
sample of weak Mg II systems span the redshift interval 
0.4 < z < 2.4. Photoionization constraints have already 
suggested that a range of ionization properties and metal- 
licities can be expected for the low ionization phase. It 
would be unusual to assume that the entire population 
of weak Mg II systems are tracing some unique type of 
physical process/structure, given these variations and the 
large redshift interval surveyed. 

To investigate, we compared the observed Fe II to 
Mg II ratio between the various systems, as it is a reliable 
constraint on density and chemical enrichment history. 
Figure [T7l shows the measured N(Fe ll)/iV(Mg ii) as a 
function of redshift. Because of the many non-restrictive 
limits at small A^(Mg n), particularly for low redshift 
clouds, it is hard to evaluate whether there is a signifi- 
cant relationship between N(Fe n)/iV(Mg n) and z. In 
order to consider this issue, we separated clouds with log 
A^(Mg ii) < 12.2, those that were likely to have only lim- 
its on Fe II (based on inspection of Figure [8]), and consid- 



ered only the stronger of the weak Mg n clouds, plotted 
in red color in the top panel of Figure [TT] It appears that 
there is an anti-correlation between N(Fe n)/A^(Mg ii) 
and z. At high redshift, there is an absence of de- 
tections with larger A^(Fe li)/7V(Mg n) values, while 
at low redshift there are many detections with large 
N(Fe n)/iV(Mg ii) and few limits that could even be 
consistent with small N(Fe n)/iV(Mg n) values. We 
applied a Kolmogorov-Smirnov (K-S) test to compare 
the distributions of N(Fe ll)/JV(Mg ii) at z > 1.5 and 
z < 1.5 for the clouds with log N(Mg ii) > 12.2. In the 
cases where only upper limits are available, we conserva- 
tively include these as values when performing the K-S 
test. The distributions are shown in histogram form in 
the lower panel of Figure[17] We find that there is a prob- 
ability of only P(KS) = 0.006 (KS statistic D = 0.505) 
that the two samples are drawn from the same distribu- 
tion. The probability is likely to decrease if upper limits 
could be replaced with actual detections. Thus we find 
that the observed anti-correlation of N(Fe n)/7V(Mg n) 
with z is statistically significant for log AT(Mg n) > 12.2 
clouds. 

We have found an absence of log N(Mg n) > 
12.2 clouds at high redshift with large values of 
N(Fe n)/7V(Mg ii) and an apparent absence of log 
AT(Mg ii) > 12.2 clouds at low redshift with small val- 
ues of N(Fe ll)/N(Mg ii). There are a number of low 
redshift clouds with limits that could be consistent with 
small values. Furthermore, for weaker clouds (with 
log iV(Mg ii) < 12.2), there are some examples of low 
N(Fe n)/JV(Mg n) values at low redshifts. We conclude 
that large AT(Fe n)/7V(Mg ii) clouds are present only at 
z < 1.5 and not at higher redshifts. The other popu- 
lation with small N(Fe ll)/JV(Mg n) exists both at low 
and high redshifts. 

As demonstrated earlier, systems in which log 
A^(Fe Il)/iV(Mg ii) > -0.3 are constrained to have a 
high density (log U < -3.7, n H > 0.05 cm -3 ). Those 
with a lower Fe II to Mg II ratio have lower densities, 
ranging down to log U < —2.0, which corresponds to 
tih < 0.001 cm -3 . The observed trend in the Fe II to 
Mg II ratio with redshift, therefore, could imply the ab- 
sence of high density clouds in the low ionization phase 
in weak absorbers at high-z. Such variations in the phase 
structure are plausible if the weak systems are probing a 
different combination of astrophysical systems/processes 
at z ~ 2 and z ~ 1. Furthermore, if the absorbers are op- 
tically thin H I clouds, then we are also seeing a change 
in the thicknesses of the low ionization gas clouds, from 
kiloparsec-scale at z ~ 2 to a range of values including 
both parsec-scale and kiloparsec-scale clouds at z ~ 1. 

Alternatively, gas clouds that are enriched primarily by 
Type II SNe events will have [a/Fe] > 0, in which case 
the observed Fe II to Mg II column density ratios will be 
low. Thus the observed trend could also indicate that 
the weak Mg II clouds are predominantly a-enhanced at 
high redshift, with an increasing contribution to the pop- 
ulation at lower redshift from clouds with a higher iron- 
abundance. Increasing the [a/Fe] in the Cloudy models, 
would then lead to low Fe II to Mg II clouds having high 
densities (log U < —3.0, tih > 0.01 cm -3 ), similar to 
the iron-rich clouds. The relevance of abundance pat- 
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tern variations is discussed in detail in § [HI where we 
speculate on the physical origin of these absorbers at the 
two redshift epochs. 

6. WEAK Mg II ABSORBERS & SATELLITES OF STRONG 
Mg II SYSTEMS 

Strong Mg II systems are understood to be absorption 
arising in the disk and extended halos of normal galaxies 
([Bergeron et al. I HMD : ISteidel et al. 111991 l20f)l . Their 
broad (Av ~ 150 km s _1 ) and kinematically complex 
Mg II li ne profiles are found to be c onsistent with this 
picture (|Charlton fc Churchill I [1998) . A characteristic 
feature in many z ~ 1 strong Mg II systems is weak, kine- 
matic subsystems separated i n velocity from the dom - 
inant absorption component ([Churchill fc Vogtl [2001). 
Such kinematic subsystems are likely to be gas clouds 
in the extended halo of the absorber in an arrangement 
analogous to the Galactic high velocity cloud (HVC) and 
intermediate velocity cloud (IVC) populations. In ear- 
lier work, we hypothesized that a non-negligible frac- 
tion of weak Mg II systems could be the extragalactic 
analogs of Milky Way HVCs, in which a random line 
of sight intercepts the surrounding halo cloud(s), but 
misses the optically-thick absorber. The possibility of 
such an event is favored strongly by some recent obser- 
vations which find a patchy distribution (less than unity 
cov ering factor) for the gas in the extended halos o f galax- 
ies (|Tripp fc Bowen Il2005t [Churchill et al. Il2007f) . For a 
patchy halo, a sight line that passes only through the 
clouds in the halo is more li kely to produce a weak Mg n 
system than a strong one ([Churchill et al. I [2005). Our 
hypothesis was primarily based on the observed evolu- 
tion in the redshift number density (dN/dz) of weak 
Mg II systems and the evolution in the gas kinematics 
of strong Mg II ab sorbers over the same redshift interval 
of 0.4 < z < 2.4 ([Narayanan et al. 1 120071 : iMshar et all 
I2007T) . 

To extend this postulate further, and also to test its 
validity, we compared the Fe II to Mg n ratio for the 
low ionization gas in weak Mg n systems to that for 
the satellite clouds of st rong Mg n systems presented in 
iChurchill fc Vogtl (l200l) . Based on an observed break in 
the velocity distribution of Voigt Pro file components in 
their sample of strong Mg II systems, IChurchill fc Vogtl 
(2001) specified clouds at |Aw| > 40 km s" 1 as interme- 
diate velocity or high velocity subsystems (i.e. satellite 
clouds). The satellite clouds in that sample were sep- 
arated in velocity by as much as \Av\ ~ 350 km s -1 
from the system center, with a median value of |Au| = 
165 km s -1 . In Figure HH we plot the Fe n to Mg n 
column density ratio of these subsystems and compare 
it to the same in our sample of weak Mg II clouds. We 
have included only those weak absorbers that are within 
0.4 < z < 1.2, equivalen t to the redshift interval of the 
IChurchill fc Vogtl ([200li) sample. 

The comparison shows that the weak Mg n clouds 
closely resemble the satellite clouds of strong Mg II sys- 
tems. To begin with, the column density of Mg n in the 
satellite clouds spans roughly the same range of values 
as that of weak Mg n systems. The Fe n to Mg n in the 
satellite clouds have a scatter which is also comparable 
to the scatter in weak absorbers at the same redshift. 



A Kolmogorov-Smirnov test estimates that the two sam- 
ples are consistent with being drawn from the same dis- 
tribution [P(KS) = 0.633, D=0.196]. Comparable to the 
subset of iron-rich weak absorbers are several satellite 
clouds with 7V(Fe n) ~ A(Mg n), which consequently 
constrains their density to njj > 0.05 cm -3 . In addition, 
a significant subset of the satellite clouds also have much 
lower densities (njj < 0.001 cm -3 ), which make them 
analogous to the a-enhanced weak Mg n clouds. 

7. SUMMARY 

Using a large sample of recently discovered weak Mg n 
systems ([Narayanan et aL~1l2007l ). we have derived con- 
straints on the chemical and ionization conditions in their 
low ionization gas. In addition to Mg II, we have mea- 
sured the equivalent widths and column densities of a 
number of other prominent metal lines associated with 
these absorbers. The significant results reported in this 
paper can be summarized as follows: 

1. In our sample of 100 weak Mg II systems, we find 
that only 48% are single cloud abs orbers. This f r action 
is smaller than the past results of iRigbv et al. I (2002) 
where the majority (67%) of weak Mg n absorbers were 
found to be single cloud systems, but is consistent within 
errors. The VLT/UVES sample that we consider in this 
paper is a fact or of ~ 5 la r ger th an the Keck/HIRES 
sample used bv lRigbv et al. I (|2002f ). We find no evidence 
for an evolution in the ratio of single to multiple cloud 
absorbers over 0.4 < z < 2.4. 

2. We find the equivalent widths and column densities 
of C II and Si n are well correlated with the equiva- 
lent widths of Mg II, with minimal scatter in the re- 
spective relationships. The column densities of C n and 
Si II yield the following relationships with Mg II; log 
N(C n) = (0.82 ± 0.10) log N(Mg n) + 3.42, and log 
N(Si n) = (1.02 ± 0.25) log A(Mg n) - 0.26. The pres- 
ence of a significant correlation in the equivalent widths, 
extends the possibility of using C n and Si n as proxy 
doublets for detecting analogs of weak Mg II systems at 
z > 2.5 in the optical spectra of quasars. 

3. If a large fraction of weak Mg n clouds are sub- 
Lyman limit systems (i.e. optically thin in H I with 
N(H i) < 10 17 cm -2 ), then the observed column density 
of Mg II constrains the metallicity in the low ionization 
gas to Z > 0.1 Zq. We also find the neutral fraction 
of magnesium to be very low in almost all weak Mg n 
clouds, approximately ~ 2 dex smaller than the corre- 
sponding 7V(Mg n). 

4. From assuming a solar abundance pattern, we find 
that the clouds for which 7V(Fe n) ~ A(Mg II) have their 
ionization parameters constrained to log U < —3.7, cor- 
responding to uh > 0.05 cm -3 . If the low ionization gas 
is optically thin in neutral hydrogen, then this places an 
upper limit of R < 10 pc on the thickness of these gas 
clouds. Similarly, clouds with A(Fe n) << A(Mg n) 
are constrained to have higher ionization parameters 
(log U ~ —2 in some cases) and lower densities. If 
the weak Mg n clouds, in which Fe n is observed to 
be weak relative to Mg II, are a-enhanced, then that 
would yield higher constraints on density similar to the 
7V(Fe n) ~ N(Mg n) absorbers. 
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5. In the past, ionization models using CLOUDY (ver- 
sion 90) have often not succeeded in recovering the ob- 
served Mg I to Mg II ratio in both strong and weak Mg n 
systems. The ionization fraction of Mg I, compared to 
Mg II, predicted by the models was not sufficiently large 
to explain the observed TV (Mg i)/7V(Mg n). Therefore, a 
separate, cold (T ~ 500 K), high density (njj > 1 cm~ 3 ) 
phase, centered at the same velocity as the Mg II phase 
was proposed in order to recover the observed Mg I in 
the models. However, in the current version of Cloudy 
(ver 07.02.01), with improvements in the rate coefficients 
of charge transfer reactions, the model Mg I to Mg n 
fraction is higher by ~ 0.5 dex for a given ionization pa- 
rameter log U . Such an increase makes it consistent for 
Mg I to be in the same low ionization phase as Mg II, in 
both weak as well as strong Mg n systems. 

6. Most of our C II and Si n measurements are for 
systems at z > 1.5. In single phase models, the con- 
straints from C II and Si n are typically high for the 
ionization parameter (log U > —2.5), which is incon- 
sistent with the constraints derived for clouds in which 
N(Fe n) ~ iV(Mg n). However, we also find an evo- 
lution in the relative strength of Fe II, compared to 
Mg II, such that towards higher redshift (z > 1.5) there 
might be a paucity of iron-rich systems (see Figure I17|) 
. The absorbers in our sample, for which there is si- 
multaneous coverage of C II, Mg II and Fe II, suggest 
that the N (Fe II) could be sufficiently small compared to 
AT(Mg n) in the high redshift clouds. Moreover, a non- 
negligible fraction of C II can arise in the high ionization 
gas, traced by C in and C IV, such that C II, in itself, 
cannot be used to determine the physical conditions in 
the low ionization gas in weak absorbers. 

7. We find that deviations from a solar abundance pat- 
tern is required to explain the observed column density 
of Al II in many weak Mg II clouds. In particular, sys- 
tems in which N(A\ n) < 0.17V(Mg Ii) require the abun- 
dance of aluminum in the low ionization gas to be lowered 
by ~ 0.7 dex, consistent with a-enhacement. Models 
with super-solar metallicity generally produce less Al II 
relative to Mg II, but some reduction of the aluminum 
abundance is still required for many clouds. When the 
abundance of aluminum is reduced, models underpredict 
Al III absorption unless the ionization parameter is high, 
which is sometimes inconsistent with that derived from 
other ions. This suggests that Al III, like C II, sometimes 
arises partly in a separate, higher ionization phase. 

8. In our sample, we find a relative absence of weak 
Mg II clouds with iV(Fe n) ~ N(Mg n) at high redshift 
(z > 1.5) compared to many detections of A^(Fe n) ~ 
AT(Mg n) towards low-z. This observed trend can be 
interpreted in two ways : (1) an absence of high den- 
sity, low ionization gas at high-z and/or (2) the pres- 
ence of [a/Fe] > in weak Mg n clouds at high- 
z. The other population of weak absorbers, in which 
N(Fe n) << N(Mg n), are detected at all intervals 
within 0.4 < z < 2.4. 

9. We find similarities between the observed column den- 
sity of Mg II as well as the Fe n to Mg n column density 
ratio in weak Mg II clouds and the high velocity sub- 
systems (i.e. satellite clouds) of strong Mg n absorbers. 
The range of N(Mg n) and N(Fe ll)/iV(Mg n) for the 



two groups are comparable. This could be suggestive of 
the fact that some fraction of weak absorbers could be 
probing a similar type of physical structure as the satel- 
lites of strong Mg n systems. 

8. DISCUSSION 

Weak Mg II absorbers have bee n identified over a large 
redshift interval < z < 2 .4 (|Churchill et al. 1 119991 : 
iNaravanan et al.~ll2005l I2007T ) . corresponding to a great 
majority of the history of the universe. Within this in- 
terval, their redshift number density (dN/dz) is found 
to be evolving, wi th a peak value of dN /dz = 1.76 
at < z >= 1.2 (jNaravanan et al. I l2007f l Towards 
lower redshift, the decrease in number density follows 
the expected curve for a non- evolving population (fo r 
a ACDM concordance model) (|Naravanan et aL~i r2005). 
At z > 1.2, the dN/dz has been found to decrease 
rapidly, such that an extrapolation to z > 3 would yield 
a value of zero. In other words, the observed redshift 
number density does not suggest that a significant pop- 
ulation of weak Mg II systems exists at z > 3. In con- 
trast, the number density of Lyman limit systems (LLSs) 
has been found to increase towards high redshift. At 
z ~ 0.7,1.5 and 3, the dN/dz o f LLSs is estimated to 
be 0.7, 1.1 and 1.9, respec tively (|Stengler-Larrea et al. I 
I1995L ISargent et al. I Il989l ). These values are in turn 
closely matched by the redshift number density of strong 
Mg II systems (with 1^(2796) > 0.3 A) at those same 
redshifts (e.g. iNestor et al. 1 I2005D . Therefore, a sub- 
stantial fraction of the observed weak Mg II clouds at 
z ~ 0.7 and z ~ 1.5 ought to be gaseous structures 
that are optically thin in H I (i.e. sub-Lyman limit sys- 
tems with iV(H i) < 10 17 cm~ 2 ). This, consequently, 
would constrain the metallicity in the low ionization 
gas of many weak Mg II absorbers to Z > Z Q , in or- 
der to reproduce the observed column density of Mg II. 
Detailed photoionization models, where information on 
the H I column de nsity has been available, further sup- 
port this inference (ICharlton et al~l l2003; Ma siero et al. I 
120051 : iMisawa et al. Il2007h . 

Bearing in mind these observed number statistics and 
constraints on the chemical and ionization conditions 
described in this work, we now proceed to discuss the 
plausible hosts of these low ionization, high metallicity 
weak absorbers. Given the range of ionization proper- 
ties and chemical abundances, it would be unusual to as- 
sume that the entire population of weak Mg II systems 
would correspond to one unique type of astrophysical 
process/structure at all redshifts. 

iSchave et al. I (|2007h have recently suggested that weak 
Mg II clouds are likely to arise in gas ejected from star- 
burst supernova-driven winds during an intermediate 
stage in free expansion, before settling in pressure equi- 
librium with the surrounding IGM. The galaxy popula- 
tions detected at high redshift (z ~ 2 — 3) are found to be 
rapidly star-formi ng, with a star-formation rate of 10 — 
100 Mq yr" 1 fe.g. lPettini eraT1l27)0il IChoi et al. Il200fih . 
The starburst events associated with these could give rise 
to galactic scale outflows that can displace large amounts 
of chemically enrich ed gas from the ambient ISM into the 
exten ded halo (e.g. Hcckm an et al. 1 12001: Pc ttini et al. I 
2001). The strong clustering of C IV systems with Ly- 
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man break galaxies, which dominate the star formation 
rate at high-z, is possibly a si gnature of such outflows 
(|Adelberger et al. 1 12003. 2005a). Such supernova driven 
winds are observed to have a multiphase structure, with 
a non-negligible fraction of the interstellar gas in a warm 
neutral phase (T < 10, 00 K) traced by such l ines as 
C II, Si II, Fe II and Al II (jSchwartz et al. II2006T) and a 
cold neutral component (T ~ 100 K) de tected in Nal 
(jHeckman et al. 1 120001: fRupke et al. 112002ft . A sight line 
that directly intercepts the outflow close to the star- 
burst region is likely to produce a very stron g, saturated, 
and k inematically broad absorpt ion feature (iBond et al. I 



2007ft . 



as 



1200 1| ). However, as described in lSchave et al 
the wind material moves farther into the outskirts of the 
extended halo of the galaxy, the column densities would 
decrease in response to a decreasing density in the am- 
bient medium. At this stage, fragments in the wind, 
generated through hydrodynamical instabilities, would 
manifest as weak Mg II clouds, and later as weaker C iv 
absorption associated with H I lines in the Ly-a forest 
(jZonak et al. 112004 ISchave et al. 112007ft . 

The interstellar clouds, ejected from correlated super- 
nova events, are likely to be highly chemically enriched 
because of the close association wit h the feedback from 
star formation. ISimcoe et al. I (|2006ft discovered evidence 
for such chemically enriched gas (Z > 0.1Z©) at z ~ 2.3, 
at distances of ~ 100 — 200 kpc from luminous star- 
forming galaxies, which they interpret as feedback from 
supernova winds or perhaps tidally stripped gas. The 
low ionization lines such as Mg n, Fe n, Al II, C II, and 
Si II in the absorbers presented in that study have col- 
umn densities similar to those of weak Mg n clouds in 
our sample. Material that is directly related to star- 
forming events is likely to have [a/Fe] > 0. Weak Mg n 
clouds associated with such events would therefore have 
N{Fe n) << iV(Mg n). This is consistent with the dom- 
inant population of high redshift (z ~ 2) weak Mg n 
clouds, and with some fraction of the clouds towards low 
(z ~ 1) redshift. The high metallicity weak C IV ab- 
sorption clouds presented in ISchave et al. I (|2007i ) were 
estimated to have sizes that are small (R ~ 100 pc), less 
than the Jean's length for self-gravitating clouds, imply- 
ing that they arc likely to be short lived. Such a transient 
physical nature is also a feature of weak Mg n clouds 
(Narayanan et al. 2005), and is anticipated for relics of 
winds. 

The metal enriched interstellar gas expelled from the 
disk would resemble the high velocity gaseous structures 
surrounding t he Milky Way , as th ey move through the 
galaxy's halo. Elli son et al. I (|2004ft , from estimating the 
coherence scales of low and high ionization gas associated 
with weak absorbers, have suggested a scenario in which 
weak Mg n absorption could arise in the outskirts of or- 
dinary galaxies, where the filling factor of the low ion- 
ization clouds (i.e. number of clouds per cubic parsec) is 
small compared to that in the center. For low ionization 
gas, the coherence scale is ~ 2 kpc, i.e., there is a high 
probability of seeing weak Mg n absorpt ion along two 
lines of sight separated by this distance (|Ellison et al. I 
l2004ft . However, the separate low ionization phases must 
not fully cover this ~ 2 kpc region, since individual ab- 
sorbing clouds are not seen along bot h lines of sight sepa - 
rated by tens to hundreds of parsecs (jRauch et al. Ill999h . 



Also, photoionization models have shown that cloud line- 
of-sight thicknesses are often < 100 pc. This suggests a 
clustering of separate clouds on a ~ 2 kpc scale, as well 
as implying a flattened geometry for the coherent struc- 
ture, consistent with the findings of iMilutinovic et al. I 
(2006) . This scale coul d be consistent with dwarf galax- 
ies ([Ellison et al.1l2004ft or with tidal streams. Sight lines 
through gas stripped in tidal interactions of galaxies can 
also produce sub-Lyman limit systems, and related weak 
Mg II absorbers. Gas that is tidally stripped in merger 
or accretion events could also form stars and provide a 
source of enriched gas clouds to the halos of high redshift 
galaxies. These would be analogous to the Milky Way 
circumgalactic gaseous streams, related to accretion of 
interstellar gas from satellite galaxies. 

In this context, we emphasize that the Milky Way 
analogs of weak Mg n absorbers are not likely to 
be the HVC complexes detected in 21cm and/or Ha 
emission, since those have -/ V QH i) > 10 18 cm" 2 
(jWakker fe van Woerder] 119971 : iPutman et all l2003ft . 
The weak absorbers must instead correspond to a pop- 
ulation of halo clouds with lower H I column densities. 
Spectroscopic observations in the ultraviolet along var- 
ious sight lines through the Milky Way halo have de- 
tected such high velocity gas in which the H I column 
densities are sub-Lyman limit (N(H i) ~ 10 16 5 cm" 2 



iCollins et alMooi iFox et al.ll2005l : iGaneulv et aT1l2T)05ft . 
These clouds, which exhibit multiple gaseous phases, 
have the column densities of C II, Si n and Fe II con- 
strained to values similar to what we find for these ions 
in our weak Mg n s ample. Using t he HST/STIS archival 
spectra of quasars, iRichter et al. I (|2008ft have identified 
a population of high velocity clouds in the Milky Way 
halo in which A^(H i) < 10 18 cm" 2 , with a few having 
N(H i) < 10 17 cm" 2 . The low ionization metal lines 
associated with these halo clouds are kinematically nar- 
row and weak, identical to the high-z weak Mg II sys- 
tems. Two of the sight lines that cover Mg n measure 
Wr(2796) < 0.2 A. These observations lend further sup- 
port to the proposition that at least some fraction of the 
weak Mg n absorption systems are likely to have their 
physical origin in gas clouds residing in the halos of high- 
z galaxies. 

The star formation per co-moving volume is known 
to b e roughly constant between z ~ 1 an d z ~ 3 
(e.g. iBouwens et~al~l 120031 : IWang et al.ll2006ft . In this 
same interval however, the number density of weak 
Mg II clouds has been found to decline from a value 
of dN/d z = 1.76 at z = 1.2 to dN/dz = 0.65 at 
z = 2.2 (|Naravanan et aL~1l2007ft . If a significant frac- 
tion of weak Mg n clouds, especially those in which 
N(Fe n) << A^(Mg n), form in supernova driven out- 
flows from star- forming galaxies, then it would seem per- 
plexing that a declining trend is observed for dN/dz from 
z ~ 1 to z ~ 2.4. We would expect dN/dz of weak 
Mg II absorbers to be not decreasing so drastically if 
they were all directly connected to interactions and out- 
flows. However, we find a spread in the physical proper- 
ties of weak Mg n absorbers, and an evolution in these 
properties from z ~ 1 to z ~ 2. We would expect that 
those weak Mg n absorbers with 7V(Fe n) << iV(Mg n) 
could be consistent with an origin in superwind conden- 
sations, since that process would lead to a-enhancement 
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and could produce high metallicities. We have found 
such absorbers both at z ~ 1 and z ~ 2, and their num- 
bers are roughly consistent with a constant dN/dz for the 
sub-population, and thus with a constant star formation 
rate over the same interval. 

We then must also consider the sub-population of weak 
Mg II absorbers with N(Fe n) - N(Mg n). These ob- 
jects have a dN/dz that peaks at z ~ 1, and they are 
apparently rare at z ~ 2. When we combine the two pop- 
ulations, the 7V(Fe n) ~ iV(Mg n) and the N(Fe n) << 
iV(Mg n) clouds, the result is a dN/dz that declines 
from z ~ 1 to z ~ 2.4, as is observed. The evolution 
and the physical properties of the N{Fe II) ~ AT(Mg n) 
absorbers provide clues to their origins. It has been 
shown that, to reproduce the observed solar composi- 
tion of iron, nucl eosynthetic yields fro m Type la events 
have to included (|Timmes et al. |[T995f) . Thus, the weak 
Mg II clouds in which the Fe n to Mg n column den- 
sity ratio approximately reflects Solar abundance (i.e. 
N(Fe ii) ~ N(Mg il), see Figure |BJ) should also be Type 
la-enriched. Enrichment by Type la supernovae requires 
that there is a ~ 1 billion year delay from the onset of 
star formation until the elements produced in the Type 
la enter the interstellar medium. What is needed to ex- 
plain these clouds is a process that peaks at z ~ 1.5 
(~ 1 Gyr before the peak of dN/dz at z = 1.2) to produce 
the stars that subsequentl y give rise to the i r on en rich- 
ment. It has been no ted in lLvnch fe Charltonl (|2007f ) and 
iMisawa et al. I (|2007[ ) that this peak closely compares to 
the peak in the global star fo r mation rate in dwa rf galax- 
ies jKau^a^etaD [200l Ha^reTSZIIlSQS)- More 
generally, for Type la SNe to contribute to enrichment, 
the absorbing structure must persist for more than 1 bil- 
lion years. This may exclude superwinds from starbursts 
and tidal debris, that are likely to be relatively short- 
lived. However, it would be consistent with Type la en- 
riched gas trapped in the potential wells of dwarf galax- 
ies, or with intergalactic star-forming structures in th e 
cosmic web (jRigbv et al. ll2002tlMiIutinovic et al. Il2006f ). 
Either of these sites could also be consistent with the 
coherence len g ths o f low ionization gas as estimated in 
lEllison et all (pOOl ). 

In the loc a l uni verse, IStocke et al. I ((2004) and 
iKeenev et al. I (|2006l) have associated weak Mg n ab- 
sorbers with unbound winds from dwarf starburst galax- 
ies. In their respective examples, a post-starburst galaxy 
is identified at impact parameters of llh~ l kpc and 
33/i _1 kpc from the line of sight. Compared to mas- 
sive star-forming galaxies, the halo escape velocities are 
smaller for dwarfs and hence, they are more efficient in 
transporting metal enriched gas clouds into intergalactic 
environments. Yet, outflows associated with starbursts 
from dwarfs are also likely to be a-enhanced, such that 
the iV(Fe n) ~ 7V(Mg n) weak absorbers would remain 
largely unexplained. Thus the clouds that have large 
N(Fe n) require a site that has been previously enriched 
by local star formation, such as the dwarf galaxies them- 
selves. 

We conclude, based on the range of physical prop- 
erties derived for the weak Mg n absorbers, that they 
arise from at least two types of astrophysical processes. 
The iV(Fe n) << N(Mg n) clouds, some of which are 



a-enhanced, are produced in processes related to the 
concerted action of massive stars, such as superwinds. 
They thus evolve in number along with the global star 
formation rate in massive galaxies, so we would expect 
roughly a constant dN/dz from z ~ 1 to z ~ 5 from this 
population. The N(Fc n) ~ iV(Mg n) clouds, which 
would be enriched by local Type la supernova, could be 
housed in relatively dense pockets within the potential 
wells of dwarf galaxies or intergalactic structures. Be- 
cause the star formation rate in dwarfs peaks at a later 
epoch (z ~ 1.5), we expect these clouds to emerge from 
z ~ 2 to z ~ 1. The two populations combined would 
lead to a gradual increase in the total weak Mg n ab- 
sorber dN/dz from z ~ 2.4 to z ~ 1. This picture would 
predict that the dN/dz of weak Mg n absorbers would 
remain constant from z ~ 2.4 up to higher redshifts, 
because the star formation rate in massive galaxies was 
constant. 

However, the redshift interval 1 < z < 3 also cor- 
responds to the epoch over which hierarchical struc- 
ture growth and mergers are most active. The frac- 
tion of interacting galaxies, and proto-galaxies with ir- 
regular luminosity profiles, is observed to be higher to- 
wards high redshift (z > 1, ~ 40%. lAbraham et al. I 
19961: Ivan den Bergh et al. 1 1 19961; IConselice et al. 112003 
Elmegreen et al. 1 [2005) . Indirect evidence of this dy- 
namical evolution of galaxies was also noticed by 
iMshar et al. I (|2007f ) in the evolution of the kinematic 
profiles of strong Mg n systems. Towards high redshift 
(z ~ 2), the kinematic profiles in a number of strong 
Mg II lines were found to be particularly complex, with 
absorption in multiple clouds that were linked with each 
other continuously in velocity. This is suggestive of on- 
going accretion events. In contrast, z ~ 0.5 strong Mg n 
systems typically have a distinct region of strong ab- 
sorption due to several blended clouds, surrounded by 
one or more weaker high velocity components (i.e. satel- 
lite clouds). This is the expected absorption signature 
for quiescent disk/halo galaxies. It was proposed by 
Narayanan et al. (2007) that the superpositions of 
the numerous halo clouds present in high redshift proto- 
galactic structures lead to a rarity of isolated weak Mg n 
absorbers during that epoch. Instead, the weak Mg n 
absorbers would be consolidated into a stronger Mg n 
absorber. The result would be a deficit of weak Mg n 
absorbers at z > 2.5. Near-IR high resolution spectro- 
scopic observations are needed to determine the dN/ dz of 
weak Mg n absorbers at z > 2.5 and to determine their 
metallicities, ionization conditions, and chemical abun- 
dances. 
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TABLE 1 

EQUIVALENT WIDTH OF METAL LINES ASSOCIATED WITH 
WEAK Mg II ABSORBERS 



QSO 


A A 


z abs 


type 


Me- IT 2796 


Me- I 2853 


Fe II 2383 


Al II 1671 


C TT 1 335 


Si it 1 260 


Al in 1855 


3c336 


3530 - 


6650 


0.702901 


s 


0.028+0.004 


< 0.010 


< 0.009 










CTQ 0298 


3520 - 


8550 


1.256069 


s 


0.057+0.004 


< 0.004 












Q 0001-2340 


3060 - 


10070 


0.452394 


m 


0.138+0.003 


< 0.019* 


0.019+0.001 
















0.685957 


s 


0.035+0.001 


< 0.002 










< 0.008 








1.651484 


s 


0.077+0.001 


0.010+0.001 


< 0.004* 


< 0.003 


0.075+0.001 


< 0.081* 


0.004+0.001 


Q 0002-4220 


3160 - 


10070 


1.446496 


s 


0.042+0.000 


0.002+0.000 










< 0.001 






1.988656 


m 


0.276+0.001 




0.016+0.003 


0.081+0.001 


0.171+0.001 


0.180+0.001 


0.045+0.002 


Q 0011+0055 


3770 - 


10000 


0.487264 


s 


0.244+0.019 


< 0.040 


















1.395635 


s 


0.186+0.004 


< 0.013 










0.043+0.008 








1.777926 


s 


0.127+0.003 


< 0.051 




0.059+0.006 








Q 0013-0029 


3060 - 


9890 


0.635069 


m 


0.205+0.014 


< 0.005 


< 0.053* 
















0.857469 


m 


0.150+0.005 


< 0.005 


0.014+0.004 








0.023+0.003 








1.146770 


m 


0.047+0.001 


< 0.009 


0.004+0.001 








< 0.034 


Q 0042-2930 


3530 - 


6800 


0.798665 


m 


0.239+0.007 


< 0.007 


0.063+0.004 
















1.091866 


s 


0.162+0.005 


< 0.012 


0.69+0.001 








< 0.011 


Q 0100+130 


3520 - 


10000 


1.758694 


m 


0.108+0.008 


< 0.013 






0.125+0.004 












2.298494 


m 


0.232+0.004 




< 0.003 




< 0.291 


0.081+0.002 




Q 0109-3518 


3060 - 


10070 


0.769646 


s 


0.044+0.001 


0.002+0.000 


< 0.001 








< 0.005* 








0.896004 


m 


0.020+0.001 


< 0.001 


< 0.003 


< 0.009 














1.182696 
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0.135+0.001 


0.009+0.001 


0.016+0.001 


< 0.020* 






< 0.030* 
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0.822606 
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0.269+0.003 
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< 0.020* 








0.910117 
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0.060+0.004 


< 0.006 


< 0.003 
















1.174224 
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< 0.007 


< 0.006 


< 0.005 






< 0.005 
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< 0.009 




< 0.005 






< 0.004 
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0.199+0.006 


< 0.015 


0.020+0.002 


0.037+0.001 


0.104+0.001 
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< 0.011 
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< 0.013 


< 0.009 


0.026+0.003 






0.021+0.003 
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< 0.013 


< 0.009 


< 0.011 
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< 0.019 
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< 0.011 
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3060 - 
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0.737248 
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0.022+0.001 


< 0.001 


< 0.002 








< 0.003 






1.708492 
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0.026+0.001 


< 0.002 
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0.009+0.001 
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0.168+0.008 


< 0.013 
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1.269042 
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< 0.016 


< 0.015 


< 0.031 






< 0.017 


Q 0329-2550 


3060 - 


10070 


0.992899 
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0.279+0.002 


< 0.001 


0.027+0.002 








0.005+0.001 








1.398230 
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0.025+0.001 


< 0.006 


< 0.003 










Q 0329-3850 


3070 - 


8500 


0.929608 
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0.073+0.002 


< 0.003 


< 0.003 


< 0.007 






< 0.004 








0.970957 
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0.055+0.001 


< 0.003 


< 0.004 


< 0.011 






< 0.031* 


Q 0429-4901 


3050 - 


10080 


0.584249 
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0.016+0.002 


< 0.004 


< 0.004 
















1.680766 
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0.015+0.001 




< 0.012 


< 0.009 


0.012+0.001 


< 0.024 


0.005+0.001 


Q 0453-4230 


3060 - 


10070 


0.895865 


s 


0.034+0.001 


< 0.001 


< 0.001 


< 0.003* 














1.039517 


m 


0.189+0.003 




0.105+0.001 










Q 0549-213 


3500 - 


6640 


1.343495 


s 


0.181+0.010 




< 0.012 


< 0.12 






< 0.013 


Q 0551-3637 


3060 - 


9370 


0.505437 


m 


0.052+0.014 


< 0.013 


















1.491972 


s 


0.176+0.004 


< 0.018 


< 0.010 


0.043+0.005 


0.108+0.010 




< 0.019 


Q 0810+2554 


3050 - 


6640 


0.821727 


m 


0.252+0.001 


0.019+0.001 










< 0.003 








0.831511 


m 


0.158+0.002 


0.004+0.001 




< 0.008 






< 0.031* 


Q 0926-0201 


3060 - 


10000 


1.096336 


s 


0.020+0.001 


< 0.001 


< 0.006 


< 0.013 






< 0.023 








1.232206 


m 


0.069+0.004 


< 0.003 


0.010+0.001 


< 0.002 






< 0.005 


Q 0940-1050 


3110 - 


10070 


2.174535 


s 


0.035+0.001 


< 0.008 


< 0.005 


0.004+0.000 


< 0.029* 


< 0.033* 


0.011+0.000 



TABLE 1 — Continued 





QSO 


A A 


z abs 


type 


Mg II 2796 


Mg I 2853 


Fe II 2383 


Al II 1671 


C II 1335 


Si II 1260 


Al m 1855 


Q 


1122-1648 


3060 - 


10070 


0.806215 


m 


0.249+0.001 


< 0.003 


0.032+0.001 








0.045+0.001 










1.234140 


m 


0.200+0.000 


0.010+0.000 


0.017+0.000 


< 0.133 






0.028+0.000 


Q 


1140+2711 


3775 - 


10000 


2.196632 


m 


0.193+0.002 


< 0.061* 






0.173+0.001 


< 0.295* 




Q 


1151+068 


3705 - 


10000 


1.153704 


s 


0.108+0.003 


< 0.008 


< 0.004 








< 0.016 


Q 


1157+014 


3520 - 


7400 


1.330502 


s 


0.120+0.002 


< 0.007 


0.020+0.004 


0.024+0.006 






0.013+0.002 


Q 


1158-1843 


3070 - 


10070 


0.506041 


s 


0.022+0.001 


< 0.003 


< 0.021* 


















0.818119 


m 


0.063+0.001 


< 0.001 


0.009+0.001 








< 0.004 


Q 


1209+0919 


3520 - 


7770 


1.264983 


s 


0.083+0.007 


< 0.019 


< 0.015 










Q 


1229-021 


3530 - 


6650 


0.700377 


s 


0.010+0.001 


< 0.004 


< 0.006 


















0.768862 


s 


0.033+0.002 


< 0.005 


< 0.005 


















0.830858 


m 


0.134+0.000 


0.011+0.001 


0.032+0.002 










Q 


1418-064 


3765 - 


9945 


1.516673 


s 


0.075+0.003 


< 0.019* 










< 0.047* 










2.174224 


s 


0.178+0.004 


< 0.023* 


0.018+0.003 










Q 


1444+014 


3520 - 


5830 


0.509719 


m 


0.193+0.007 


0.022+0.003 


< 0.129* 


















1.101989 


m 


0.033+0.005 


< 0.007 


0.005+0.000 








0.032+0.005 


Q 


1448-232 


3060 - 


10070 


1.019191 


m 


0.033+0.005 




< 0.002 


< 0.070 






< 0.002 










1.473201 


m 


0.269+0.009 


0.011+0.000 


0.007+0.001 


0.024+0.001 




< 0.116* 


< 0.051 










1.585464 


m 


0.088+0.002 


0.008+0.001 


0.005+0.000 


0.005+0.001 


0.069+0.001 


0.044+0.001 


0.022+0.002 


Q 


1621-0042 


3530 - 


6800 


1.174521 


m 


0.237+0.012 


0.004+0.001 


0.003+0.001 










Q 


1629+120 


3050 - 


6800 


1.379330 


m 


0.142+0.007 


< 0.018 


< 0.017 










Q 


2000-330 


3495 - 


9945 


1.249864 


s 


0.032+0.001 


< 0.003 


< 0.004 










Q 


2044-168 


3520 - 


9900 


1.342525 


m 


0.057+0.004 


< 0.013 




< 0.014 






< 0.013 


Q 


2059-360 


3750 - 


9280 


1.242973 


s 


0.015+0.001 


< 0.006 










< 0.009 










1.399947 


m 


0.109+0.002 


< 0.006 


0.006+0.002 








< 0.011 


Q 


2116-358 


3530 - 


6640 


0.539154 


s 


0.115+0.004 


< 0.018 




















0.775358 


m 


0.183+0.010 




< 0.013 










Q 


2132-433 


3500 - 


6640 


0.793600 


m 


0.184+0.007 


0.012+0.005 


0.030+0.007 










Q 


2204-408 


3520 - 


6800 


1.335279 


m 


0.052+0.004 


< 0.009 


< 0.008 










Q 


2206-199 


3420 - 


6640 


0.948363 


m 


0.265+0.007 


< 0.005 


0.016+0.002 








< 0.007 










1.297044 


s 


0.148+0.001 


0.008+0.001 


0.005+0.001 


< 0.064* 






< 0.005 


Q 


2217-2818 


3060 - 


9890 


0.599512 


s 


0.115+0.003 


0.005+0.001 


< 0.016* 


















0.786515 


m 


0.204+0.001 


< 0.001 


0.025+0.005 








< 0.028* 










1.054310 


m 


0.046+0.002 


< 0.002 


< 0.002 


< 0.021 
















1.082780 


m 


0.125+0.001 


< 0.001 


< 0.001 


< 0.003 






< 0.003* 










1.200162 


m 


0.099+0.002 


< 0.001* 


< 0.002* 








< 0.022* 










1.555849 


m 


0.268+0.001 


< 0.076 


0.043+0.000 


0.036+0.000 


< 0.271* 


0.125+0.002 


< 0.006* 


Q 


2222-3939 


3530 - 


6640 


1.227553 


s 


0.114+0.005 


< 0.019 


< 0.013 


< 0.013 






< 0.012 


Q 


2225-2258 


3050 - 


10000 


0.831374 


m 


0.033+0.002 


< 0.003 


< 0.005 


< 0.03 






< 0.005 










1.433018 


m 


0.167+0.002 


< 0.011 




0.036+0.003 


0.162+0.004 


0.010+0.010 


0.009+0.002 


Q 


2243-6031 


3140 - 


10000 


0.828081 


m 


0.263+0.005 


< 0.003 


< 0.007 








< 0.009 








1.389707 


m 


0.106+0.022 


< 0.003 




< 0.007 
















1.755704 


s 


0.108+0.001 


0.007+0.001 


< 0.004 


< 0.008 


< 0.106* 


< 0.058* 




Q 


2314-409 


3520 - 


6640 


0.843114 


s 


0.043+0.003 


< 0.008 


< 0.010 










Q 


2347-4342 


3100 - 


10070 


1.109640 


s 


0.040+0.004 


< 0.003 


< 0.003 
















1.405367 


s 


0.074+0.001 


< 0.003 


< 0.004 




< 0.095* 




< 0.030* 










1.796233 


s 


0.160+0.001 


< 0.009 


0.010+0.001 




< 0.163* 


0.066+0.001 


0.004+0.000 



Note. — * indicates lines that were contaminated by absorption features at other rcdshifts. In most cases, the contamination was from H I lines of the Ly-a forest. The second column gives 
the wavelength coverage of the UVES spectrum of each quasar, the third column gives the redshift of each absorber, and the fourth column indicates whether the weak Mg II absorption was in a 
single cloud (s) or in multiple clouds (m). Columns 5-11 arc the total rest-frame equivalent widths of the respective lines. 
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NOTE: Table 2 lists the velocity, column density and Doppler parameter for Mg I, Mg II, Fe II, C II, 
Si II, Al II and Al ill associated with the 100 weak Mg II systems in our sample. The table (18 pages) 
can be downloaded from http://www.astro.wisc.edu/~anand/vptable.pdf 
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Fig. 1. — A few example systems plots from the 100 weak Mg II absorbers in our sample represented on a velocity scale with the system redshift corresponding to km s _1 . The 
vertical axis is continuum normalized scale for the spectra. The vertical tick marks correspond to the centroid of the absorption components in each spectral feature, derived from Voigt 
profile fitting. Regions that are blueward of the Ly-a emission of the quasar are labeled in the respective panels as in Ly-a forest. The quasar name and redshift of the system are also 
labeled in each system plot. (Note: The full set of system plots will be presented in the online version of the journal.) 
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Fig. 2. — The top panel shows the frequency distribution of the number of weak Mg II clouds (i.e. number of Voigt profile components) 
in our sample. Single cloud systems account for 48% of the total number of systems. The bottom panel compares the rest-frame integrated 
equivalent width of Mg IIA2796 line as a function of the number of clouds. Most of the weaker systems [W r (2796) < 0.1 A] have absorption 
only in one or two clouds. 
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Fig. 3. — Distribution of single and multiple cloud weak Mg II systems as a function of redshift. The redshift distribution is comparable 
between the two types of systems (N c = 1 and N c > 2; i.e., single and multiple cloud) and they do not exhibit preference for either low 
(z ~ 1) or high (z ~ 2) redshift. 
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Fig. 4. — The total (i.e. integrated) rest-frame equivalent width of Mg IIA2796 A in each absorber, as a function of the absorber's 
redshift. The red and black data points correspond to single and multiple cloud systems respectively. It is evident that there is no 
significant correlation between the strength of weak absorber and its redshift. 
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Fig. 5. — The comparison of the total rest-frame equivalent widths of Mg n (WV(2796)) in each absorber, with other corresponding low ionization transitions and Al III. The large solid 
dots are 3<x detections, the large open circles are detections that are affected by blending with an absorption feature at some other redshift, in which case the measurement is considered 
as an upper limit, indicated by a downward pointing arrow from the open circle. Non-detections at the 3cr level are plotted using just the downward pointing arrow. In most cases, the 
blending was from H I lines of the Ly-a forest. The equivalent width corresponding to non-detections at the 3<r level are plotted using downward pointing arrows. The dash- dot line in 
the C II and Si II panels are the y = x line for comparison with Mg II. The correlation coefficients (Spearman's p) between VF r (Mg n) and IVrCFe II), W r (Mg i), W r (C II), Wr(Si II), 
WV(A1 Ii), WV(A1 III) are 0.64(0.00), 0.50(0.00), 0.75(0.00), 0.38(0.20), 0.67(0.00), and 0.38(0.00) respectively. The value in parentheses represents the significance level, i.e., the probability 
that the observed value of p would be greater than or equal to the actual value by chance. 
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Fig. 6. — The comparison of Mg II column densities with the column density of other low ionization transitions and Al III. The large solid dots are 3a detections, the large open circles 
are detections that are affected by blending with an absorption feature at some other redshift, in which case the measurement is considered as an upper limit, indicated by a downward 
pointing arrow from the open circle. The column density of non-detections, estimated from the 3it equivalent width limit, are plotted using just the downward pointing arrows. The solid 
line is a linear regression fit formalizing the relationship between the two ions. The standard deviation of the corresponding fits are indicated by the dashed-dot line. The regression lines 
are dra wn for the three ions (Fe II , C II and Si II) that display t he lea st sc atter with N(Mg n). T he red dashed line in each panel indicates the solar abundance pattern based on values 
given in Grcvcssc & Sauval (1998); Allcndc Pricto ct al. (2001, 2002) and Holwcgcr ct al. (2001). In the lower left panel, the red dashed line therefore corresponds to y = x. 
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Fig. 7. — Cloudy photoionization curves indicating how the column density of Si II, C II, Mg II, Fe II, Al II and Mg I change with ionization parameter (log U), metallicity (Z) and 
neutral hydrogen column density [iV(H i)]. The density (njj) was calculated using the expression log njj = log n 7 - log U, for log n 7 = -4.70, corresponding to the number density of 
ionizing photons [hv > 13.6 eV) at z = 2. 
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Fig. 8. — Cloudy grid of photoionization models for O.IZq, Zq and IOZq metallicity with measurements of Fe II and Mg II column density over-plotted. Vertical curves correspond 
to lines of constant JV(H i) and the horizontal curves lines of constant ionization parameter (log U). The Cloudy models were computed for the intensity of the extragalactic ionizing 
background radiation (EBR) at z = 2 and z = 1. Weak Mg II clouds at z > 1.5 are plotted in the EBR at z = 2 panel and clouds at z < 1.5 are plotted in the EBR at z = 1 panel. The 
large solid dots are 3cr detections, the large open circles are detections that are affected by blending with an absorption feature at some other redshift, in which case the measurement 
is considered as an upper limit, indicated by a downward pointing arrow from the open circle. The column density of non-detections, estimated from the 3cr equivalent width limit, are 
plotted using just the downward pointing arrows. 
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Fig. 9. — The distribution of Fe II to Mg II column density ratio for single and multiple cloud systems. The histogram filled with the light 
color includes measurements that are upper limits in Fe II, whereas the histogram filled with the dark color excludes those. For multiple 
clouds, we have included the Fe II to Mg II ratio in individual clouds. It is evident from the figure that both single and multiple cloud 
systems span roughly the same range of Fe II to Mg II column density ratios. 
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Fig. 10. — Cloudy grid of photoionization models with measurements of Mg I and Mg II column density over-plotted. The description of the Cloudy curves and the data points are the 
same as in Figure [S] 
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Fig. 11. — Cloudy grid of photoionization models with measurements of C II and Mg II column density over-plotted. The description of the Cloudy curves and the data points are the 
same as in Figure[5] Almost all weak Mg II absorbers with coverage of C IIA1335 are at z > 1.5 and therefore we plot all of them on the z = 2 EBR plot. 




Fig. 12. — Cloudy grid of photoionization models with measurements of Si II and Mg II column density over-plotted. The description of the Cloudy curves and the data points are the 
same as in Figure[S] Almost all weak Mg II clouds with coverage of Si IIA1260 are at z > 1.5 and therefore we plot all of them on the z = 2 EBR plot. 
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Fig. 13. — Cloudy grid of photoionization models with measurements of Al II and Mg II column density over-plotted. The description of the Cloudy curves and the data points are the 
same as in Figure [3] 
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Fig. 14. — Cloudy grid of photoionization models with measurements of Al III and Mg II column density over-plotted. The description of the Cloudy curves and the data points are the 
same as in Figure [3] 
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Fig. 15. — Left Panel : The figure compares the C II to Mg II ratio with Fe II to Mg II in weak Mg II clouds that have measurements for both C II and Fe II. The filled circles correspond 
to weak Mg II clouds that had firm measurements for both C II and Fe II. The open circles correspond to data points in which either C II or Fe II or both have measurements that are 
upper limits. Middle Panel : compares the Al III to Mg II ratio with Al II to Mg II in weak Mg II clouds that have measurements for both Al III and Al II. Right Panel : compares the 
C II to Mg II ratio with Si II to Mg II in weak Mg II clouds that have simultaneous measurements for both C II and Si II. 
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Fig. 16. — On the left panel is the Al III to Al II ratio in weak Mg II clouds discussed in this paper. Measurements that are upper limits in Al in and Al II are indicated using downward 
and upward pointing arrows respectively. The right panel shows the Al III to Al II ratio in sub-DLA and DLA systems taken from the literature. If weak Mg II clouds are optically thin 
in H I (i.e. iV(H i) < 10 17 cm -2 ), then the measurements indicate that the anti-correlation between 7V(H i) and Al III to Al II ratio that is observed for DLA and sub-DLA systems, also 
extends to lower N(H i) values. 
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Fig. 17. — Top panel : shows the Fe II to Mg II ratio in all the weak Mg II clouds as a function of redshift of the system. The 
filled circles are firm detections (i.e. detections that are not limits), the large open circles with arrow pointing downward corresponds to 
Fe II lines that are affected by blending with an absorption feature at some other redshift, in which case the measurement is considered 
as an upper limit. The non- detections at the 3cr level are plotted using downward pointing arrows. Those weak Mg II clouds in which 
N(Mg II) > 10 12 ' 2 cm -2 are plotted in red. Bottom panel : shows the distribution of Fe II to Mg II in absorbers at z > 1.5 and z < 1.5. 
The width of each bin in the distribution is 0.3. The data used to create the frequency distribution includes only those weak Mg II clouds 
in which iV(Mg n) > 10 12 ' 2 cm -2 , below which a large fraction of systems only have upper limits for Fe II. 
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Fig. 18. — The figure compares the column density ratio between Fe II and Mg II for weak Mg II clouds and the high velocity components 
(i.e. satellite clouds) of strong Mg II systems. In the right panel the downward pointing arrows indicate measuremen ts that are upper 
limits in 7V(Fe II). The Fe II and Mg II column densities for the satellite clouds of strong Mg II systems is taken from Churchill &: Vogti 
( 2001), and are represented as large blue points. The small red points correspond to weak Mg II clouds presented in this paper. The 
left panel compares the distribution of the column density ratio between Fe II and Mg II for weak Mg II clouds and the satellite clouds 
of strong Mg II systems. The blue histogram is for the satellite clouds that are offset in velocity from the main absorbing component by 
|Av| > 40 km s" 1 . The histogram shaded with lines is for the satellite clouds that are offset in velocity by [Au| > 90 km s _1 . The red 
histogram is the distribution for the sample of weak Mg II clouds presented in this paper. 



